I — Introduction 



The Moyal phase space representation provides a self-contained and complete formula- 
tion of quantum mechanics, wherein the statement of dynamical evolution closely parallels 
that of classical Hamiltonian mechanics. In the Moyal formalism^ each quantum observ- 
able A is represented, via the Wigner transform, by a unique real phase space function 
(symbol) A (q,p) and, conversely, each symbol defines an operator on Hilbert space via 
Weyl quantization. In this paper we analyze the Wigner- Weyl image of the Heisenberg 
picture evolution operator and through it obtain a semiclassical asymptotics of quantum 
evolution that is caustic free. These non-perturbative expansions are valid to arbitrary 
order in Planck's constant, are applicable to physically interesting observables, and hold 
for systems having a Hamiltonian symbol that is a smooth function on phase space. 

Let A(t, s) be a quantum observable in the Heisenberg picture, for a system having an 
arbitrary time dependent Hamiltonian H(t). If a : A \— > A w = <r (A) denotes the linear 
Wigner transform^ 2-8 ! (from operators in Hilbert space H = L 2 (R d ) to functions on the 
phase space T*R d ) then the evolution equation obeyed by the Weyl symbol A (t,r) is 

| : A w (*,r) = -{A w (*,r),^ w (r)} M . (1.1) 

Here the Moyal bracket { , } M is the image of the quantum commutator times (iK) -1 under 
the Wigner- Weyl isomorphism. When applied to (1.1) Bohr's correspondence principle is 
realized in a transparent fashion. The Moyal equation of motion (1.1) has the classical 
(% = 0) limit 

■^A(t,r) = -{A(t,T),H c (r)}, (1.2) 

wherein { , } is the Poisson bracket. The fact that (1.2) is a transport equation for the 
classical dynamics governed by Hamiltonian function H makes it easy to build h J, 
expansions for quantum observables using the classical flow in phase space. Moreover, it 
turns out that equations (1.1), and (1.2), are naturally accompanied by initial conditions 
on A , and A; this feature is responsible for the caustic-free nature of the attendant semi- 
classical expansion. We continue this introduction by providing a more detailed overview 
of this and other results in the paper. 

In the usual Hilbert space operator based description of quantum mechanics, the uni- 
tary evolution U(t, s) generated by H(t) is the operator valued solution of the Schrodinger 
initial value problem 

ih^-U(t, s) = H(t)U(t, s), U(s, s) = I; (1.3) 
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s is an arbitrary initial time, and / is the identity in 7i. A Schrodinger picture observable 
A is transformed to the Heisenberg picture by the linear mapping r (from operators to 
operators) which conjugates A with U, 

r(s,t) : A i-> U(t,s)i AU(t,s) = A(t,s). (1.4) 

The Schrodinger operator U(t, s) evolves states, whereas its Heisenberg counterpart -T(s, t) 
evolves observables. Either of these two operators carries all dynamical information about 
the system and each is equally basic. However, the Heisenberg picture evolution operator 
r(s,t) has a much simpler semiclassical behavior than does U(t,s). For this reason we 
focus our study largely on r(s,t). Good reasons for writing the time labels of r in the 
order (s,t) will become apparent later. 

Consider the ingredients that enter the Wigner-Weyl representation of r. 

As noted above, each observable A is identified with a function A w on phase space 
through the Wigner transform a . (Subsection 2a reviews the Wigner-Weyl method.) An 
observable is called semiclassically admissible if its Weyl symbol is asymptotically regular 
at h = 0, 

oo . r 

^(0 = ^(0 + ^^(0. (i.5) 

r=l 

Here the variable ( = (g,p) denotes a generic point in the 2<i-dimensional phase space 
R 2d ~ T*R d . The leading term of the series (1.5), A c ((), is regarded as the classical 
counterpart of A. It is easy to check that the common observables of quantum mechanics 
such as position, momentum, angular momentum, kinetic energy, the Hamiltonian and 
algebraic combinations of these operators are semiclassically admissible. If an observable 
depends on additional variables, such as the times (t,s), then these labels are appended 
to all the functions in (1.5), e.g. A w becomes A (t,s). 

The classical analog of U(t, s) is the measure preserving phase flow g(t, s \ ■ ) : R 2d — > R 2d 
which solves Hamilton's equation 

±g(t,s\0 = JVH c (t,g(t,s\0) (1-6) 

with initial condition g(s,s\() = ( £ R 2d . In (1.6) J is the Poisson matrix carrying the 
symplectic structure of Hamilton's equation. 

When the Heisenberg evolution r(s,t) is viewed in the Wigner-Weyl representation, 
it becomes the composition V s t = a w o r(s,t) o cr" 1 . Stated less abstractly, an equivalent 
characterization of T s t is the following: Suppose A w (t,r) is the solution of (1.1) with the 
t = t initial value A ; then 

T s y.A w ^AJt,s) (1.7a) 
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is the counterpart of (1.4) expressed as a linear mapping among phase space functions. 
The action of T s t implements time evolution for the Weyl symbols of Heisenberg picture 
observables. We will write this mapping in the functional form 

r M K)=A w (M), (1.76) 

and refer to T s t as the Heisenberg- Weyl evolution operator. A key feature of T s t is that 
this mapping has a regular behavior with respect to Planck's constant h = 2ttU } and admits 
a small h formal asymptotic series 



n=0 



which is of a semiclassical nature. The leading term is the pull-back (on functions) of 
the classical phase flow, 7^ = g(t,s \ ■ )* . Thus for an arbitrary time dependent function 
a(t, ■ ) on phase space, A(t, r, • ) = j^a(t, • ) = a(t, g(t, t\ ■ )) solves the transport equation 
(1.2). 

Two different methods of computing the expansion (1.8) will be presented. The first 
method, found in §3, determines the evolved Weyl symbol of an arbitrary operator A(t) by 
first expressing it as a function of the canonical phase space coordinate operator z = (q,p). 
The Heisenberg equation obeyed by Z(t,s) = r(s,t)z, 

j t Z(t,s) = JVH w {t,Z(t,s)), (1.9) 

is a quantum version of Hamilton's equation (1.6). In (1.9) the necessary operator ordering 
prescription is understood to be Weyl ordering [see (3.6)]. The operator solution of (1.9) 
preserves the canonical commutation relations for all time; in Weyl symbol form: 

{z*jt,s),zi(t,s)} M = rJ. (1.10) 

It is evident from (1.9) and (1.10) that it is appropriate to interpret Z(t,s) and its Weyl 
image Z (t,s) as the definition of a quantum trajectory. Furthermore, it will be shown 
that this quantum trajectory carries a generic semiclassical expansion whose leading term 
is the classical flow g. 

With the aid of a new cluster-graph representation of the symbol of an exponenti- 
ated operator, it is possible to efficiently solve (1.9) for the the h expansion coefficients 
of Z (t,s). The higher order quantum fluctuations about the classical flow obey inho- 
mogeneous Jacobi (variational) equations which may be solved recursively using a phase 
space Green function adapted to initial conditions. These basic ingredients are then used 
to construct the h expansion of A w (t } s), from which it is straightforward to identify series 
(1.8). 
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Given the striking and simple similarities between the quantum and classical equations 
of motion (1.1) and (1.2), it is worth recalling some of the difficulties that have restricted 
the wider use of the Moyal formalism. The equation of motion (1.1) has an unpleasant 
non-local form. Specifically, when written in terms of phase space derivatives the bracket 
{ , } M contains derivatives of arbitrarily high order. Thus (1.1) is a pseudodifferential 
equation in 2d + 1 variables. Similarly, the Heisenberg equation (1.9) for the canonical 
variables Z(t, s) is a 2d system of coupled ODE's which (for arbitrary Hamiltonians) is 
infinite-order non-linear with solutions that are non-commuting unbounded Hilbert space 
operators. These are substantial difficulties. A good portion of this paper is devoted 
to developing new calculational methods and representations that are effective in this 
environment. 

The most novel of these representations is the above-mentioned cluster expansion, 
which provides a potent addition to the Wigner-Weyl calculus. Derived in §2b, the cluster 
expansion has the form 

oo 1 

(^) w (0 = exp{^-£ J (C^)}, (l.H) 

where the cluster functions C- are determined by summation over simple connected graphs 
on j vertices. Each vertex represents the symbol A w and the graph edges are Poisson 
bracket operators coupling these symbols. Representation (1.11) has wide applicability. 
The case where i is a static (time independent) Hamiltonian and e A is its evolution 
operator is discussed in §2d. The central role played by the exponential operator e w ' z in the 
Wigner-Weyl method makes this a natural arena for the use of (1.11), as is demonstrated 
in §3. 

A second method for computing T s t is presented in §4. There, its Moyal equation 

§; T T,t = -{rT,t,H w (T)} M (i.i2) 

is solved directly, order-by-order in h, after (1.8) is substituted as an ansatz. The resulting 
recursive formulae for 7^ are very compact and do not require previous computation of 
the Z (t,s) expansion. A comparison of the first few semiclassical expansion coefficients, 
generated by the two different approaches, shows that they are mutually consistent. 

This second approach has antecedents^ 9-11 ^ which consider the evolution of quantum 
densities (rather than observables) for static Hamiltonians. The difference between the 
Heisenberg equations for density and observable evolution is merely a sign. However, 
when time dependent observables and Hamiltonians are considered, as they are here, it 
turns out, somewhat unexpectedly, that the appropriate evolution problem results from 
the 'backward' rather than the 'forward' equation (in (1.1), d/dr instead of <9/<9t), as is 
discussed in §4a. 
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There are two basic limitations imposed on the type of physical system considered 
here. First, it is assumed that the wave functions are scalar valued and consequently that 
the relevant quantum state space is L 2 (R d ). Second, we presume that the classical phase 
space is flat R 2d . Our results apply to systems characterized by smooth, semiclassically 
admissible Hamiltonian and observable symbols, both of which may be time dependent. 
Expansion (1.8) is constructible for t in the maximal interval about s for which the flow 
g(t, s \ ■ ) is defined on R 2d . For smooth Hamiltonians with controlled (*-growth, the flow is 
globally defined for all time. 

Although in principle one might obtain all our results by staying strictly within a 
framework intrinsic to the Weyl symbol calculus, we have found it useful to employ both 
the operator calculus and the symbol calculus, treating both points of view on an equal 
footing. The mathematical style of most of our arguments is heuristic and formal. It is 
designed to build, primarily by constructive means, the explicit semiclassical expansions 
found in Sections 3 and 4. A second generation approach to this subject would likely 
rewrite these results in the language of the Weyl pseudodifferential operators and Fourier 
integral operators. However in a first treatment these extra layers of mathematical rigor 
would tend to obscure the basic structure of the simple semiclassical asymptotics we find. 

The Heisenberg-Weyl semiclassical expansion (1.8) has a number of distinctive and 
attractive features. The operator valued coefficients 7^ are composed of the exact classical 
flow and phase space derivatives of finite order which act on the Weyl symbol of an arbitrary 
observable. The specific differential form of the coefficient is determined directly from the 
first n + 1 cluster functions C- without the need to solve any recurrence relations. If the 

(n) 

Hamiltonian of the system is Weyl quantized the expansion has even ^-parity with all 7^ / 
for odd n vanishing. 

In §5 we discuss the computation of expectation values in the Heisenberg-Weyl for- 
malism. The general features of this approach are illustrated by considering some simple 
examples. Appendix A reviews briefly the covariance of Weyl symbols under affine canon- 
ical transformations, and studies how this property is reflected in the exponential graph 
representations found in §3. Appendix B presents the connection between T s t and its 
known integral kernel representations. 

II — Weyl Symbol Calculus 

The formulation of the Weyl symbol calculus is reviewed in the first part of this section 
and stated in a manner suitable for the analysis of semiclassical dynamics. This summary 
also serves to fix our notations for quantum and classical phase space quantities. One 
cornerstone of this calculus is the *-product, which is the symbol image of the operator 
product. In §2b we develop a connected graph representation that determines the symbol 
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for the exponential of operators, namely <7 w (exp(A)). This cluster-graph expansion is 
a natural extension of the Groenewold^ 12 ] asymptotic series for the *-product. Fourier 
transform techniques are employed in §2c to extend the cluster representation to describe 
<7 W (/(A)) for suitably smooth functions /. As a sample application of these combinatorial 
graph identities, the Weyl symbol of the Schrodinger evolution operator exp(—itH/U) is 
discussed in §2d together with the appropriate realization of the Hamilton-Jacobi equation 
which serves as the basis of a WKB expansion for this propagator. 

2a. The Wigner—Weyl Correspondence 

In classical Hamiltonian mechanics, the state space of a system with configuration 
manifold R d is the cotangent bundle T*R d . The Cartesian chart q on R d induces canonical 
coordinates ( = (q,p) = (q 1 , ■ ■ ■ ,q d ,Pi , • • -Pj) on this 2<i-dimensional phase space, which 
is then regarded as R 2d ~ R^xR^. The geometry of phase space^ 13 ! is determined by 
the symplectic form u = dp i A dq l . Of greater import here will be the Poisson structure 
determined by the matrix J inverse to that of u. The Poisson bracket of two phase space 
functions / and g is the function {/,<?} defined by 



{f,g}(0 = v/(CWV0(C), J 



o 8 
-8 



(2A) 



where ( now represents a generic point in R 2d , and 8 is the unit matrix (here (/-dimensional). 

Our discussion of the Weyl calculus in this section is not intended to be complete, 
regarding both the list of topics covered and the level of mathematical rigor. For example 
we do not discuss the statistical physics motivations which lead to the Wigner transform (cf. 
§4c and Refs. [1,2 ,7,12]), and we assume throughout that various formulae have meaning 
when the operators and symbols appearing in them are chosen from appropriate spaces. 

For a variety of reasons, it is useful to associate to each operator A (in some class) a 
corresponding phase space function (perhaps generalized) called its symbol, and denoted 
cr(A). There are many different ways^ 6 ' 9 ' 14 ' 15 ! to set up such a symbol <->■ operator cor- 
respondence. Most of these share a basic collection of common features. For example, 
it is reasonable to require a linear bijective correspondence, and that operator functions 
of q alone, or p alone, have as their symbols the same function of q } or p. The symbol 
of the 2d- vector operator z = (q,p) will then be (*, understood as a linear polynomial. 
Acting on scalar wave functions tp £ L 2 (R d ), the d position operators q = (q 1 , . . . , q d ) are 
defined by q l ij>(x) = x l ij)(x), and the conjugate momentum operators p = (p 1 , . . . ,p d ) by 
Pjij)(x) = —iUdtp I dx J . 

Consider a less trivial function, r((*), that is a polynomial of degree two or greater. 
Then, generally, there are several different operator polynomials in z which reduce to r(() 
if z is formally replaced by (. They differ in their ordering of the non-commuting operators 
q\ p- . For a particular definition of the correspondence, only one of these operators can 
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have r as its symbol. This factor ordering aspect of the symbol <->■ operator correspondence 
is an unavoidable consequence of the canonical commutation relations 



[z k , z l ] = ih.J kl . (2.2) 

The correspondence induced by normal ordering of q and p is familiar from classical 
PDE theory. There, if the operator A can be presented in the standard form 



j3£W d 
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where f3 is a d- multi-index and W the set of non-negative integers, then it is given the 
symbol a = a s 

j3£W d 

which is called the 'standard' or 'g-before-p'^ 9 ' 16 ] symbol. 

The Wigner-Weyl correspondence results from a more symmetrical choice of order- 
ing of the components of z, which we now summarize. A general phase space function 
A(() is 'quantized' (roughly, ( is replaced by z) by defining 'A(i)' according to the Weyl 
prescription:^ 17 ] in the Fourier resolution 

A(C) = [ dua(u)e iu ' c (2.3a) 

it is meaningful to replace ( by z, and the resulting operator 

A = [ dua{u)e tw - z (2.36) 



may be regarded as A(z). The Wigner-Weyl correspondence takes A <->■ A(z) as an 
associated symbol <->■ operator pair. In particular, note that e*™'^ <-> e in ' z . 

More precisely, the Weyl symbol of operator A, denoted by <r (A) or simply A w , 
is defined as the Fourier transform, (2.3a), of the unique function (or distribution) a for 
which representation (2.36) holds. One should note at this stage that in general the symbol 
A w can possess a parametric dependence on h. 

Conversely, a given function A w will be the Weyl symbol of the operator 



A = (2vr)— J du J d(AJ()e^^>, (2.4) 

obtained by inverting (2.3a) and substituting into (2.36). Formula (2.4) relates A to A w 
via integration. In formal calculations it is often very convenient to employ a differential 
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form of this relation. This is gotten by writing A w (()e 
and performing the (^-integral. The result is 



1, 



M = 



[K^uY~ lu< ] ^ (2.4) 



(2.5) 



Formula (2.5) clarifies the type of symmetrical ordering involved in Weyl quantization. 
Specifically, the basic symbol A w (() = ( a (a G W 2d ) is seen to be associated with 
the 'symmetric product' of operators^ 11 ] (|a| !) _1 V^(u • . For example, one has the 
association^ 18 ] 



(i'}"i,P,r 



Of course, an operator can be characterized by functions other than phase space-based 
symbols. A prime example is its integral kernel, i.e. the Dirac matrix element ( x | A | y ), 
for which the following formulae are useful. Taking the (x \ ■ ■ ■ \y) matrix element of (2.4) 
leads to a construction of the kernel starting from the symbol, 



(x\A\y) =h~ d / dpe ip - {x - y) l % A w ( 



x + y 



,P 



(2.6) 



(in deriving this one needs ( x | e^™ 1 ' M2 )' z | y ) = e tWl ' ( - x+y ^ 2 S(x — y + Tiu 2 ).^ One passes in 
the opposite direction, from kernel to symbol, via the Wigner transform^ 



A w (q,p) = J ^e-W^ + I, 



A 



q v 

H 2 



(2.7) 



This familiar formula, often taken as the definition of the Weyl symbol, results from (2.6) 
by changing variables and inverting the Fourier transform. It follows readily from (2.7) 
that a (A^) = a (A); in particular self-adjoint operators have real symbols. 

If the symbol representation of quantum mechanics is to be advantageous, it is neces- 
sary to have a complete and practical symbol calculus. Of special relevance is the symbol of 
an operator product. For this basic structure there are again both integral- and derivative- 
based formulae which are useful in varying circumstances; for z E U 2d they are 



(AB)Jz) = (nhy 



-2d 



dCdCAJz + OBjz + C)ef<- J <' : 



(AB)Jz) = exp { -V z • JV z ,} A w (z)B w (z')\ z ,__ 



(2.8) 
(2.9) 



The two formulae (2.8) and (2.9) have a rather different mathematical status. The 
former constitutes a composition rule constructing the Weyl symbol of operator AB in 
terms of the symbols of A and B. Under rather general hypotheses^ 3 ' 19 ~ 22 ^ on the operators 
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A and B this formula will be an exact, rigorous identity. In particular, (2.8) admits non- 
smooth symbols if the integration is understood in a distributional sense. The second 
statement — Groenewold's formula (2.9) — is a consequence of (2.8) and provides a small 
h expansion of (AB) . Hormander^ 20 ] and others^ 23-27 ! have attached various rigorous 
asymptotic interpretations to the expansion (2.9) by restricting both symbols A w , B to 
suitable spaces of smooth functions. For computational purposes, it is the algebraic and 
analytic form of the Groenewold formula that is its most important feature. 

The right side of identity (2.8) may also be viewed as defining a bilinear, non- 
commutative, associative product of symbols. Denoted by * and written 

A,*B w = (AB) w , (2.10) 

this product is sometimes referred to as the starj 28 ' 29 ! twisted^ 8 ' 25 ' 26 ! or Weyl^ 16 ] product. 
An elementary, but important, consequence of expansion (2.9) is that the *-product of two 
semiclassically admissible symbols is again semiclassically admissible. 

Differences of operator products, i.e. commutators, are fundamental in quantum me- 
chanics. Their image in the Wigner-Weyl formalism is the Moyal bracket { , } M defined 
by 

, Bj M = A w * B w - B w * A w = [A , B] w . (2.11) 

Just like the commutator, the Moyal bracket is bilinear, skew and obeys the Jacobi identity. 
The Moyal bracket appearing in (1.1), for example, corresponds to the commutator in the 
Heisenberg equation of motion. Clearly, both the *-product and Moyal bracket operations 
on phase space functions are defined with reference to a specific value of the parameter h. 

To verify formulae (2.8) and (2.9) efficiently, begin by computing AB using represen- 
tations of the form (2.36). After using the Baker-Campbell-Hausdorff (bch) formula^ 30 ] 
to combine the exponential operators one has 

^-/*,/*,«(„)K W )e-*.-.«-.-, (2,2) 

where B is the Fourier transform of 6, as in (2.3a). Note that the quantity ihJ has arisen 
from the canonical commutator (2.2). Now, (2.12) expresses AB in a form essentially 
parallel to (2.36), so the relation analogous to (2.3a) is, immediately, 

(AB) (z) = j du j dwa(u)b(w)e-^ u - Jw e l(u+w) - z . (2.13) 



To obtain (2.8), substitute into (2.13) the Fourier representations of a, 6 in terms of A w , B 
to find 



(A6)Jz) = fd(f dC A JOB JC) f j^e 



D iu-(z-C) e iw(z-C') e -fu-Jw_ 
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After shifting the ( and (' integrations by z, one can do the u-integral to obtain S(( + 
hJw/2). This delta function trivializes the remaining ^-integral and (2.8) results after 
minor rearrangement. 

Next, to obtain (2.9) from (2.13), simply re-write the exponential factors in the inte- 
grand as 



exp{jV z .JV ') 



1 i(u - z+w z') 



In the circumstances where the exponentiated derivative may be taken outside the u- and 
^-integrals of (2.13), these integrals reduce to A (z) and B (z 1 ), respectively [cf. (2.3a)]. 
This gives (2.9). It is interesting to note how the 7i-coupled derivatives in (2.9) originate 
directly from quantum non-commutativity; if q and p were to commute then the symbol of 
AB would be the ordinary product A w B w . It is the interchange of limiting order, related 
to pulling the differential operator V z • <7V ; through the du dw integration, that turns the 
result into an asymptotic expansion rather than an identity. 

In many real applications it is not enough to consider the product of just two operators. 
Arbitrary products of operators, and functions of an operator, arise. Formula (2.9) is a 
basic first step for handling such cases. While the ensuing subsections will develop a more 
complete calculus, it is helpful to introduce here the concepts and notations which allow 
generalization of (2.9) to an arbitrary product. 

Thus consider n operators A^ (/£n = {l,...,n}) with corresponding Weyl symbols 
Aff . We introduce a product function defined on n copies of phase space, 

n 

XA^ : (R 2d ) n ^ C, X_A«(C 1 , , C J = U A^(Q. (2.14) 

In the case where all the symbols are the same, we often write simply X — A w . If n is small 
it is convenient to use a 'tuple' notation for the product function and write, for example, 
a( 2 K = X -AW 

' w ^ — l£2 w ' 

Let V • be the phase space gradient striking the j th argument of this function, 

v, x_a£)](c 1 ,...,C„)= [vaW]^)]!^ ^). 

Notice that the V -'s mutually commute. In association with the symbol product X A^ 
we introduce an extended Poisson bracket operator B- , which acts on the variables (*■ and 
(■ and is defined by 

In this terminology the standard Poisson bracket (2.1) for a pair of functions / and g 
is recovered as {/, <?}(C) = B 12 ^f ■, 9^(0 where the single argument ((*) signifies the 
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diagonal evaluation , ) | _^ _ . Restated with this new notation the Groenewold 
formula (2.9) reads 

(^ (1) ^ (2) ) w (C) = e fBl2 ^i 1) ,4 2) KC). (2-16) 



Lemma 1. The Weyl symbol of a product of n > 2 operators has the form 




(2.17) 



Proof: If in (2.16) is replaced by A^A^ then 

(A (D A ( 2 ) A (3) )w(C) = e fVr^ A (i) (()(A ( 2 ) A (3) )w(0 |^^ 

= e fVc^V c , e fB 23A a )(C)Al2)(C , )A (3) (C , ) |^^ 

Now V^, reduces to V 2 + V 3 and (2.17) follows for the case n = 3, and for general n by 
iterating this type of argument. ■ 

The connection between the Moyal bracket, defined in (2.11), and the classical Poisson 
bracket is easily seen using (2.16). Since J is skew, B 21 = —B 12 , and so for a pair of smooth 
phase space functions A and H one computes 

i%{A, H} M = e^ Bl2 ^A, Hy - e^ Bl2 <H , Ay = (e^ Bl2 - e^ B21 )^A, Hy, 
K H} M = f sin(^ 12 HA, Hy = ± { ~l^ B^^A, Hy. (2.18) 

The semiclassical form of the Moyal bracket operator is thus { , } M = { , } + 0(Ti 2 ). This 
formal Poisson bracket limit is used in passing from (1.1) to (1.2). The symbol of the 
anti-commutator AH + HA results from reversing the relative sign in (2.18), to give 
z cos (f B 12 )^A, Hy. The resulting cos +z sin decomposition of the *-product fornW 31 ] 
the basis of a superalgebra. 

Another beneficial feature of the Wigner-Weyl calculus is its natural interaction with 
phase space translations, or more generally with affine canonical transformations. These 
properties are reviewed in Appendix A. 
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2b. Cluster Expansion of (e 



w 



Because of its central role in many calculations, it is necessary to have a general 
formula that determines the Weyl symbol of the exponential of an operator. We use a 
graph-combinatorial method to achieve this goal. 

The following Lemma will be useful below. Notationally, Yl^rn denotes summation 
over all multi-indices m = (m 1 , . . . , rn n ) £ W n which partition the positive integer n, i.e. 
obey the constraint 

n 

Y^jm j = n. (2.19) 

We also define = 1. 



Lemma 2. (Exponential Faa di Bruno) Let <f> be a complex valued function of a real or 
complex variable, which is C n in the neighborhood of point x. Then 

/ 7 \ n n 

1-^) e+W =e^Y, (n)c r n Il( dJ ^ x )) mj ( 2 - 20 «) 

where 

c m = n\ [ mi \m 2 \ ■ ■ ■ m n !(l!) mi • • • (n!) m "j _1 . (2.206) 



Proof: This follows from the well-known Faa di Bruno formula for the n derivative 
of a composite function. See Ref. [32]. ■ 

The following graph terminology is required by the cluster expansion theorem below. 
The set C- denotes all simple connected graplW 33 ! (clusters) which may be formed on the 
vertex set j = {1, 2, . . . , j}. Thus C £ C- means C = {j, E} } where the edge set E consists 
of unordered pairs a of distinct integers in j, linking a connected graph. A graph is simple 
if its edge set has at most one edge (or 'link') between each pair of vertex elements. A 
graph is connected if the edge set contains a unbroken path of links between each pair of 
vertex elements. The two integers comprising the edge a will be denoted Aa = mina and 
Va = max a. 

Theorem 1. Let A be an operator with Weyl symbol A w . Then the symbol of e A has 
the exponential cluster expansion 

oo _ 

( e ^ w (C) = exp{^-£ j (C^)} (2.21a) 
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where C- is the cluster function 

V = E [ II ( e(in/2)Ba - !)] >< A w(0, (2-216) 

and the Poisson bracket operator for link a is defined by 

B a = V Aa -JV Va . (2.21c) 

Proof: Using formula (2.17) with each = A yields 

oo _. oo _. 

(«' i ).=E^"). = Esi II H + %]><-4 w (2.22,.) 

n=0 n=0 l<i<j<n 

where 

L {j = e (tn/2)Bi i - 1. (2.226) 

Definition (2.226) systematically splits the effect of the Groenewold operator exp[(iU/2)B i j] 
into its pure classical, h independent, piece (the identity) and the remaining non-classical 
contribution which is consolidated into L--. Since the operators B-- commute with one 
another so do the operators . The product J^f^ <J - [1 + L-] is a polynomial in whose 
terms may be identified with simple graphs on n = {1, . . . , n}. Such a graph has edge 
{z, j} if and only if link operator appears in the term. In general this graph is a union 
of clusters (possibly trivial) on subsets of n. If I C n, define the cluster sum operator 

s i= e h l « ( 2 - 23 ) 

CeC(I) a£E 

where L a = L Aa y . Let V = {I x } represent a partition of n into disjoint subsets I x . 
Then 

IT 1 + hi\ = E II S Ix (2.24) 

the sum being taken over all possible partitions V. 

When (2.24) is applied to X — A w , as required by (2.22a), then for each fixed V the 
result factors according to the subsets I x , 

ii[ i +^]>< a w=e n { s i,* a ^- ( 2 - 25 ) 

i<j n v i x ev Ix 

The value of each factor here depends only upon the cardinality of I x . Thus consider the 
standard j-element (1 < j < n) subset J and define the j th order cluster function in accord 
with (2.216), 

Cj = S J XA w . (2.26) 

j 
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Then (2.25) can be written 



i<j n T j=l 



where m - {V) denotes the number of subsets I x in V of cardinality j . Clearly m ■ = m - {V) 
obeys (2.19). 

Next, notice that the summand in (2.27) only depends on V via the 'occupation 
numbers' \m-{V)\. Various distinct partitions possess the same multi-index m{V) = m. 
For a given ra, the number of such partitions is found to be simply c m , defined in (2.206). 
Combining this fact in (2.27) with (2.22a) gives 

oo _ n 
n=0 m j=l 

The n = term here is 1, as required by (2.22a). 
Define the power series 

OO A 

, . . \ ^ X J „ 



J- 



which reduces to the argument of the exponential in (2.21a) when x = 1. Clearly C- 
d J <f>(0) and so application of the Faa di Bruno formula (2.20a) in (2.28) gives 



(A 



"-ii" (^>) 



x = 



n=0 

by the Taylor series representation of . ■ 

The exponential cluster representation given in Theorem 1 is a simple and powerful 
computational tool, and is one of the main results of this paper. The combinatorial aspects 
of the proof of Theorem 1 are identical to those of the cluster expansions that are known for 
the statistical mechanical grand partition function, I 34 ! and for the quantum propagators of 
perturbed quadratic Hamiltonians J 35 ^ The algebraic framework of (2.21) is notably distinct 
from those applications of the cluster method, because it has no attendant integration for 
each symbol, and the basic binary coupling structure is skew rather than symmetric. 

Two properties of the cluster functions C- should be pointed out now. First, it is clear 
that C- is complex j-linear in the Weyl symbol A w . Second, notice that the link operator 
L a = e ( lfl / 2 ) B a> — 1 is formally 0{h r ) } and that in (2.216) each cluster C (being connected) 
has at least j — 1 links a. Thus 

£ J ((;n) = 0(V- 1 ) (2.29) 

provided that A w = O{Ti ) (as is the case for semiclassically admissible observables A). 
The first three cluster functions are particularly simple, 

h) = A w ((), £ 2 (C; h) = L 12 -<A w , A w y, (2.30) 

£ 3 (C; %) = {(L 12 L 23 + L 13 L 23 + L 12 L 13 ) + (L 12 L 13 L 23 )}^A w , A w , A w y((). 
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2c. The Symbol of f (A) 

The results in §2b are tailored to exponential functions, but their scope can be broad- 
ened by introducing a generalized Fourier transform. Consider functions representable in 
the form 



f(x) = J d(i(a) e cax , (2.31) 

where /i is a real or complex measure on R, and c £ C. If (2.31) has meaning for x in the 
spectrum of A, one can define in a natural way the operator 



f(A) = j d(i(a)e caA . (2.32) 

We now proceed formally and take the Wigner transform of (2.32) and employ Theo- 
rem 1, 

f(A)jC)= I d^a)ex P {f^ ( -^ J C j (C;n)}. (2.33) 
J j=i 3 ' 

The explicit extraction of (ca) J here means that C- is the cluster function for A as defined 
by (2.216). Since C 1 = A w , (2.33) may be manipulated into 



/°° jC f) °° jC 

<W«) exp { £ -f( -)■} eH = ex-p { £ -± 8*}/(A.) 

j=2 •> ' y i/ w ^._ 2 J • 



where 9 is ordinary differentiation of /. A final useful rearrangement is to perform a 
cumulant expansion of the exponential, to obtain a series in higher derivatives of /. This 
is facilitated by Lemma 2 (with x replaced by 9), and the result is 

f(A)JO = f(AJQ) + E E II Tfr^TT K'-f(AJO). (2.34) 

n=2 ™, \ ; = 2 V - J j- I 

mi=0 V / 

This shows that the Weyl symbol of operator f(A) is, to a first approximation, just the 
composition f o A w . The correction terms given by the ^2 n>2 series turn out to be 0{Ti 2 ). 
To see this, recall that C- = 0(U J X ) and further notice that the order Ti 1 contribution to 
C 2 [cf. (2.30)] is the bracket (iU/2){A w , A }, which vanishes. 
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2d. Schrodinger Evolution for Static Hamiltonians 

As a first application of Theorem 1, it is natural to consider the time evolution of 
a quantum system governed by a static Hamiltonian H . The Schrodinger evolution op- 
erator in this case is U(t) = exp(— itH /fi). Its Weyl symbol is immediately expressed in 
exponential form, 

oo 

tf w (f) = exp (2.35a) 

where C- is determined by (2.216) with A w = (t/iU)H w . From this it is seen that C- oc t\ 
and a cumulant expansion of (2.35a) will give the small t Taylor series of U (t). Of course 
the exposed dependence of A w means U (t) carries an essential singularity at ft = 0. 

If one is interested in the WKB (h J, 0) asymptotic form (rather than t — > 0) then it is 
useful to expand the link operator thus, 

i — i a 

There will be a sum over a link integer / for each edge a, so introduce the extended graph 
summation notation 

5>£(nE) 

and let r = ^2 aeE l a (r = if E = ). Then the general cluster function formula is 

£ i (*,C;n) = *^2-W _i II Tl B a a XH w (0- (2-356) 

Gj a£E a ' J 

In particular for j = 1,2: 

£l - JK H - ' ^ ~ ^ 720T 512 ' H - y - 

In the second formula, all the odd terms vanished by skewness of the Poisson matrix. 

It is also possible to investigate the Weyl symbol U (t) using WKB methods, in a 
fashion similar to the case of the familiar propagator ( x | U(t) | y ). The correct WKB 
ansatz turns out to be (at least for small t) 

Ujt, C) = exp {(z^)- x $(t, (; %)} (2.36a) 

oo 

$(f,C;ft) = £(ift)*$ fc (f,C). (2-366) 

A;=0 
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Since U(0) = /, thus U (0) = 1, and so one demands the initial conditions 

$,(0,C) = (k>0). (2.37) 

Briefly, the WKB method (which we won't carry out here) proceeds by substituting ansatz 
(2.36) into the Wigner transform of Schrodinger's equation, separating powers of h to find 
a recursive system of PDE's for {$,}, and solving them subject to (2.37) via a transport 
method. 

The leading coefficient $ (t,(*) controls the essential singularity in (2.36). It is the 
phase space analog^ 36-39 ! of the action S(t } x } y) (Hamilton's principal function), and has 
been christened the 'phase action' by Marinov in his study J 4 °] The phase space Hamilton- 
Jacobi equation which $ is required to obey via the WKB analysis is 

^(tX) = H c ((+ 1 -JV<S> (tX)). (2.38) 

Now, formulae (2.35) give a direct cluster-based construction of U (t). Let us briefly 
compare its Ti-structure with WKB ansatz (2.36). Consistency demands that (z^) _1 $ = 
Y2j>i A cluster C £ C- has at least j — 1 links and so the link integer sum r > j ' — 1. 

Thus, if H w = H c + 0(h r ), (2.356) shows that each C J = 0(n~ r ), consistent with the 
structure of (2.36). 

A next stage of comparison involves conjecturing a graphical formula for the phase 
action $ by extracting the (z^) _1 parts of each C- . To do this, one must in (2.356) replace 
H by H and keep only terms with r = j — 1. These correspond to minimally connected 
(tree) graphs, with each link integer 1=1. The resulting series for the phase action is 
remarkably simple, 



C) = tH c (() + 27^7T EI1 B « X H d0- (2-39) 

-'■ C'ETj aEE 3 



Here T- is the set of all tree graphs on j, and the j = 2 term is again absent by skewness. 

Note that the right side of (2.39) vanishes at t = as required by (2.37). Although 
it is beyond our scope to present the details here, one can also verify that it satisfies the 
anticipated Hamilton-Jacobi equation, (2.38). Since that PDE is typically infinite-order 



nonlinear in $ , this is a nontrivial result. 



Exponential cluster expansions of U (t,s) for systems with non-static Hamiltonians 
of the form p 2 /2m + v(x, r) have been derived by Dyson series-based methods in Ref. [41]. 
Euclidean as well as toroidal configuration spaces were considered there. 
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Ill — 



Semiclassical Evolution of Observables 



The first of two methods for computing semiclassical expansions for observable evo- 
lution in the Heisenberg picture and Weyl representation is presented here. We consider 
first the quantum trajectory Z(t,s). For this fundamental test case, the cluster expansion 
of Theorem 1 is used to derive Jacobi equations for the h expansion coefficients. These 
are solved in §3b using phase space Green function techniques. In the final subsection the 
expansion is extended to arbitrary observables. 

Each specific quantum system is defined by a Hamiltonian operator H(t) which gener- 
ates unitary evolution. Our derivations apply to any system described by a semiclassically 
admissible Hamiltonian; specifically, the Weyl symbol of H(t) is assumed to have the 
regular h J, asymptotic expansion 

oo . r 

H w (t,n;0 = H c (t,0 + J2^K^0. (3.1) 

r=l 

Since H(t) is Hermitian H , H c and h are all real. 

A simple but non-trivial example is that of a mass m charged particle moving in an 
external electromagnetic field, 

H(t) = ^g» v (p^ ~ a^t, q)) {p v - a v (t, q)) + </>(t, q). 

Here g^ v is a symmetric metric tensor and (a, cf>) denote the vector and scalar potentials. 
Using (2.9) to compute the symbol of the quadratic portion of H(t) readily shows that 
H (t) is just the classical Hamiltonian H (t, q } p) obtained by replacing q,p with q,p in the 
above expression. Thus no h (t, () terms appear and H(t) is just the Weyl quantization of 
H (t, (). This particular Hamiltonian has a quadratic momentum dependence. In general, 
the Hamiltonian symbols in (3.1) allow for any smooth dependence in p. 

The higher order h terms in (3.1) can originate via two inherently distinct mechanisms. 
The first one reflects the 'ordering ambiguity' in quantizing a classical observable. Different 
choices for ordering the non-commuting operators q\p- lead to different Hermitian opera- 
tors consistent with a given classical H (t, q } p) } i.e. they all have H as the U° part of their 
Weyl symbol. But the higher h coefficients differ because the commutation relations (2.2) 
permit one such operator to be expressed as another, plus discrepancy terms proportional 
to a factor of ih for each commutation. Thus, investigation of quantization schemes other 
than Weyl's will lead to Hamiltonians whose Weyl symbols contain nontrivial h r . 

The second mechanism is more intrinsic and concerns but a single operator: the true 
quantum Hamiltonian of the physical system of interest. There is no reason a priori why 
its Weyl symbol could not depend on h. This need have nothing to do with ordering issues. 
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For example, (1 + aU)p 2 /2m might be the free particle Hamiltonian, and it is ultimately 
up to experiment to detect such h dependence, if present. 

Not all common observables are Weyl quantized. For example, the Hamiltonian for 
the quantum rigid rotator has the form H = aL 2 , for some constant a and L = q X p. The 
associated symbol is H w = a(q X p) 2 — (3a/2)U 2 . 

3a. Semiclassical Quantum Trajectories 

The phase space coordinate operators z = (q } p) form a basic set of observables. In the 
Heisenberg picture they become the time dependent operators Z(t,s) = r(s,t)z, which as 
a consequence of (1.3) satisfy the Heisenberg equation of motion 

d ~ 

i% — Z(t,s) = [Z(t,s), H(t,s)] (3.2a) 



with initial condition 



Z(s,s)=z, (3.26) 



where the Hamiltonian in Heisenberg picture is H(t,s) = r(s,t)H(t). 

For each fixed s, (3.2) constitutes an operator valued initial value problem which 
determines Z(t, s), if one supposes H(t, s) is given. We will seek a semiclassical description 
of Z(t, s) by making an ansatz for its Weyl symbol 



oo . 

n 



Z w (t } s } n;0 = Y^-z r (t } s;(). (3.3) 



r=0 



Here the h dependence of the symbol Z (t,s) is explicitly displayed. The U-free coefficient 
functions z (t, s; () in the h J, formal asymptotic expansion (3.3) will be computed using 
techniques from classical mechanics. The initial conditions obeyed by these coefficients are 
a consequence of (3.26). Specifically, Z (s, s, = ( and (3.3) imply 

z (5,5;C) = C, z r ( 5 ,5;C) = (r>0). (3.4) 

The Heisenberg equation in the form (3.2a) is, in practice, not well suited for solution. 
Note that its right side contains, in addition to Z(t,s), the Heisenberg operator H(t,s) 
which is unknown. One means of circumventing this difficulty, which we shall employ here, 
is to start from the following variant of (3.2a), 

j t Z(t,s) = ^r(s,t)[z,H(t)]. (3.5) 

Here the Schrodinger commutator [i , H(t)^ is given a priori and can be expressed as a 
known function of z. Application of -T(s, t) will then produce the same function of Z(t, s), 
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thus rendering the Heisenberg equation in a suitable form. An alternative approach, used 
in §4, is to employ the 'backward' (d/ds) evolution equation instead of (3.2a). 

The symbol of the commutator in (3.5) is readily computed because the symbol 
z w (() = ( supports only one (*-derivative. Using (2.11) and (2.18) 

W , H(t)] w = ih{z> w , Hjt)} M = ihB 12 *z> w , Hjt)y = ihJ»V k H w (t). 

Now the commutator in (3.5) is represented in terms of this symbol via (2.5), 

p(t,s) = JVH w (t^ 1 -V u )e- z ^ (3.6) 

OX l u = 

This Heisenberg equation of motion was presented in the more stream-lined form (1.9) in 
order to emphasize its structural similarity to the classical Hamilton equation of motion. 
To proceed with (3.6) we swap the roles of iu and (l/z)V M using the identity 



1, 



1. 



Z(-VJ^) =g(-V u )f(u) 

l u = I 



u = 



which follows from Taylor expansion of the entire functions / and g. This gives 

|z(t,3) = e ^)- v JVif w (t,^;0), (3.7) 

where all gradients V are understood to act on the (phase space) argument of H (t, U; ■ ), 
with subsequent evaluation at ( = 0. 

Finally, take the Wigner transform of the operator identity (3.7). (in doing so note 
that V is not a Hilbert space operator (such as p), as it merely acts on the symbol H 
and can be treated like a fixed C 2d -vector for purposes of the Wigner transform.) The 
cluster expansion, Theorem 1, is applicable if we choose the exponentiated operator to be 
A = Z(t,s) • V. Thus combining (3.7) with (2.21) leads to 

d °° 1 

-Zjt, 5, %■ C) = exp{]T -£,((; %, t, s, V)} JVHjt, %■ 0). (3.8a) 

Here the dependence of the cluster functions C- upon t, s and V is shown explicitly. Thus 
the C- act as operators on the symbol H w and are determined by 

r/C; n, t, s, v) = ]T J] [ e _ x] >< z ^ a> n . () . v _ {3M) 

CeCj a£E(C) 3 

The basic objective now is to substitute the h-sevies ansatz (3.3) for Z w into (3.8), and 
then to equate coefficients of common powers of h to obtain solvable relations for {z r }^_ . 
In (3.8), ^-dependence not arising from Z w resides in the link operators L a = e ( tfl / 2 ) B a _ ^ 
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and possibly in the symbol H . Note that the gradients in the operator B act only on 
the Z functions and not H . 

W W 

Consider the classical limit of (3.8) — i.e. set h = 0. Then the link operators L a 
vanish and so C- = 0, j ' > 2, because there is no contribution from any cluster having 
links. As for C 1 , equation (2.30) shows it to reduce to 

/^(U^V) = Z w (M,0;C)-V = z (M;C)-V. 

The classical limit of (3.8) is therefore 

|z (M;C) = e Zo(M;C) - v JV# c (t,o) = JV# c (t,z (M;0). (3.9) 

Equation (3.9) is precisely Hamilton's equation (1.6) for the classical mechanical system 
with Hamiltonian H . Its solution z Q as determined by the initial condition z Q (s,s; • ) = 
Id R2 d [cf. (3.4)] is therefore the classical phase flow, g(t } s\(). That is, 

z (M; ■)=g(t,s\-):R 2d ^R 2d , (3.10) 

and, to leading order in h the quantum trajectory Z (t,s) is just the classical flow. 

In order to derive the higher h corrections z r from (3.8) it is useful to isolate this 
classical (/-based underpinning, as follows. Let the deviation of Z w from its classical limit 
g be 

6Z w = Z w -g = Y f -z r (t,s;C), (3.11) 

r=l 

which is formally 0{Ti r ). Correspondingly, C 1 splits into two parts, 

Z^M^V) =g(t,s\0-V + 8Z„(t,s,n;0-V, 
the former of which will be used to shift the argument of H in (3.8), with the result 

-Zjt, 3, h; C) = e ^w(*,-,ft;C) - v exp | J- h, t, s, V)),JVH w (*, ft; g(t,s\()) . (3.12) 

Equation (3.12) serves as a 'master equation' for extracting equations of motion for 
z r . This is done by collecting common powers of h. The h structure on its left side is 
trivially obtained via (3.3). On the right side, each of the three factors is h dependent, 
with e fi ^ w ' v and JVif w being relatively simple (cf. (3.11) and (3.1)). 

Computing z r requires analysis of C- up to order 7i . Recall from (2.29) that £j(Ci ^) = 
0(U J 1 y, this remains true for £j(Ci ^M? 3 ?^ 7 ) as given by (3.86). The equations of motion 
for the first two quantum corrections, z 1 and z 2 , will be extracted from (3.12). These 
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examples also serve to expose the form of equation arising for general r. The Ti part of 
C 2 , and the Ti 2 parts of C 2 and C 3 will therefore be needed for our calculation. 

@ 

Figure 1. Cluster C for C 2 

Consider first C 2 . In (3.86) the only cluster in C 2 is C = {2, {1, 2} }, which is illustrated 
in Fig. 1. So 

C 2 = [e^ B ^-l]^Z w -V,Z w -Vy 

= (*{zi ,zt} + ^(§)X 2 ^ , z>) v,.v fc + o(n% 

The first term vanishes because the j, k-skew matrix formed by Poisson bracket {Z^ , Z^} 
is contracted against the symmetric entity V -V fe . Inserting expansion (3.3) for Z w yields 

A = " Y B ^ ' ^ V i V * + 0(h 3 ). (3.13) 
In particular C 2 has no Ti 1 part, and so in (3.12) the sum C-/ ' j\ = 0{Ti 2 ). 

© © © 

© © © 

© © © 

Figure 2. Clusters for C 

Turn now to the Ti 2 part of cluster operator C 3 . The graphs contributing to C 3 , 
shown in Fig. 2, are three distinct labeled trees and the complete graph on 3. Since each 
L a is O^ 1 ), only the trees contribute to C 3 at the lowest (second) order. Calculation of 
this term based on (3.86) is straightforward; to within a sign each tree contributes equally 
and one finds that 

C 3 = ^B 12 B 23 ^ , gi , sHv^ + <D(h 3 ). (3.14) 
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With these results we can now extract the Ti 1 - and 7i 2 -parts of (3.12). For r = 1, 

^ Zl (M;0 = (z 1 (t,s;()-V)JVH c (t,g) + JVh^g). 
A natural and more suggestive rearrangement of this equation is 

j(t;s,C) Zl (M;C) = JV/^^MIO), (3.15) 

where J is the Jacobi operator along g, 

J(t;s,() = j t ~ JVVH c (t,g(t,s\C))- (3.16) 

(Wi? denotes the Hessian matrix of H .) This Jacobi operator is associated with vari- 
ational equations obtained by differentiating Hamilton's equation with respect to a pa- 
rameter, such as (. Thus (3.15) shows z 1 to obey an inhomogeneous Jacobi equation 
(the inhomogeneous term JV/i 1 in this instance can actually be present only for non-Weyl 
quantized Hamiltonians). 

Turning to the r = 2 case, one finds from (3.12)-(3.14) that z 2 also obeys such a 
Jacobi equation, 

J(t; s, ()z 2 = [( Zl • V) 2 - ^B 2 12 Xg-V + ^B 12 B 23 X r v] JVH c (t, g) 

+ 2(z 1 •V)JVh 1 (t,g) + JVh 2 (t,g). (3.17) 

In this instance, the inhomogeneous term need not vanish in the Weyl quantized case. 

The general features of equations (3.15) and (3.17) turn out to be generic for r > 1. 
To summarize these features, z r obeys an inhomogeneous Jacobi equation 

J(t;sX)z r = f r (t,s;(). (3.18) 

The Jacobi operator J in (3.18) arises by considering the only two terms in the h r part 
of (3.12) that contain z r : d/dt comes from the left side and JVVif c originates in the 
linear SZ w • V term in the expansion of e SZw '^. The inhomogeneous term / is built 
from g, z 1 , . . . , z r _ 1 and from H , h 1 , . . . , h . The next subsection will prove solvability 
of (3.18) when / is known. Since h i are given from the quantum Hamiltonian, this implies 
the possibility of solving the r-indexed family (3.18) recursively for {z r }. 
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3h. Jacohi Green Function Solutions 

In this subsection we explicitly construct the solution of the inhomogeneous Jacobi 
equation (3.18) subject to the relevant homogeneous initial condition (3.4), i.e. 

z r (5,5;C) = 0, r>l. (3.19) 

Equation (3.18) is a 2<i-dimensional first order ODE, with parameter ( £ R 2d . Basic 
ODE theory implies that the initial value problem (3.18)-(3.19) has a unique solution z r , 
provided the given coefficient functions / and W H (t, g) are sufficiently smooth in t. 

This problem can be solved by constructing the phase space Green function tailored 
to it. The Green function G for (3.18)-(3.19) is a matrix function G : R X R -> R 2dx2d 
obeying the fundamental inhomogeneous Jacobi equation 

J(t; s, ()G(t, p) = 6(t -p)J (pe R), (3.20) 

and the homogeneous initial condition 

G(s,p) = (p £ R). (3.21) 

The necessary dependence of G upon the parameters (s,C) appearing in J is being no- 
tationally suppressed. Once the Green function is known one verifies directly that the 
solution to (3.18)-(3.19) is 

z r (M;C)= / dpGit.pyj^f^s-Q. (3.22) 
Jr 

The construction of G based on phase space Jacobi fields is relatively straightforward. 
Recall that a Jacobi field, B } is a matrix solution of the homogeneous Jacobi equation; 
namely, B : R -> R 2dXm satisfies 

J(t;s,()B(t) = 0. 
Thus each of B } s m columns is a vector Jacobi field. 

Lemma 3. Let B be any 2dx2d matrix Jacobi field which is invertible at some time 
p . Then B(p) is invertible for all p £ R and 

G{t, p) = [Q(t -p)- G(s - p)] Bit)B(p)- 1 J (3.23) 

is an initial value Green function obeying (3.20)-(3.21). 
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Proof: A basic property of Jacobi fields is symplectic conservation: if A and B are, 
respectively, 2d X m and 2dxm' Jacobi fields, with A T the transpose of A, then 



A{t) 1 JB{t) = const. G R mXm . 

This is readily verified by ^-differentiation. Choose A = B here and take the determinant 
to find ( det B(t)) is a constant, which is nonzero by considering t = p . Hence B remains 
invertible. Next, it is immediate that (3.23) obeys (3.21), and (3.20) follows by a simple 
calculation. ■ 

To complete our construction of G a suitable Jacobi field is needed. A natural choice 
is to use the derivative of the phase flow 

B(t) = B(t; s, C) = Vg(t,s\C) = |r(MlC). (3.24) 

(The standard derivative-matrix notation is used here: dg/d( is the 2d square matrix 
with entry dg l /d( J .) That Wg is a Jacobi field for j7(t;s,C) follows by taking the 

(*-derivative of Hamilton's equation for g. Moreover at time p = s, *V g(s } s\() = 6 is 
invertible. We also note that this choice of B(t; s, () is symplectic (i.e. B T J B = J), hence 

B(p)- 1 = J- 1 B(p) T J. (3.25) 

To summarize, the solution of (3.18)-(3.19) gotten by combining (3.22)-(3.25) is 

z r (t, s; () = J dp Vg(t, s\() JVg(p, s\() T J~ V^p, s; (). (3.26) 

It should be emphasized that there are no restrictions on allowed values of t, s here, 
apart from the requirement that the classical flow be defined. This is a result of the fact 
that initial value problems are being solved. This in turn is not just a consequence of the 
initial value nature of Schrodinger's equation, but also reflects the fact that phase space 
symbols of Heisenberg picture operators are being studied. If one alternatively studies, for 
example, the conventional propagator ( x | U(t, s) | y ) or the Weyl symbol U (t,s), then 
as is well known the relevant trajectories and Green functions obey boundary conditions 
involving separate times t and s (in place of initial condition (3.21)). Such two-time 
boundary problems do not always have solutions; specifically, this occurs at conjugate 
points. The description of evolution being considered here is, by contrast, free of such 
caustics. 

We now apply the general formula (3.26) to solve in closed form the Jacobi equations 
for z 1 and z 2 derived in the previous subsection. The right side of (3.15) defines / x (t, s; () 
and so 

Zl (t, s; C) = J dp Vg(t, s\C) JVg(p, sKfVh, (p, g(p,s\()) . (3.27) 
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This simplifies by using the chain rule to express Vg T V/i 1 (/?, g) = V^/i 1 (/?,(/(/?, , and 
recognizing the resulting integrand of (3.27) to be a Poisson bracket. Thus 



Zi(M;C)= / dp{g{t,s\-),h 1 (p,g{p,s\-))}{( > ). 

J s 



(3.28) 



For z 2 we shall restrict attention to the usual case of a Weyl quantized Hamiltonian, 
H = H c . In this case h 1 , h 2 and z 1 are all zero, and combining (3.17) with (3.26) gives 

/"* r 1 

z 2 (t, s; C) = / dp Vg(t,s\()JVg(p,s\() T [ - -B 12 2 ^g a (p, s) , g\p, s)y(()V a V^ 

•J S 

+ ^B 12 B 23 ^g a ( P} s) , g P(p, s) , g\p, s)^(C)V a V /J V 7 ] VH c (p, g(p,s\0) ■ 

(3.29) 



3c. Evolution of Semiclassically Admissible Ohservahles 

The prior results of this section have succeeded in developing a recursive scheme for 
computing a semiclassical expansion of the quantum trajectory Z(t,s). Here we extend 
those results to the case of more arbitrary observables A(t, s) in a simple way. 

Let A(t) be a time dependent Schrodinger picture operator which is semiclassically 
admissible. Then its symbol has the form, fully parallel to (3.1), 

oo yr 

AJt, h; C) = A c (t, + J2^T °r(*> 0- (3-30) 

r=l 

The ^-expansion coefficients A c and a r are real if A(t) is self- adjoint. As the notation in 
(3.30) suggests, A (t, () is the classical, h independent, part of A w (t } U; (). 

In the Heisenberg picture, the evolution of A(t) is A(t,s) = r(s,t)A(t). We will 
construct the Weyl symbol of A(t, s) as an expansion of the form 

oo ^ r 

A w (M,ft;0 = ]L^a r (M;C) (3-31) 

r=0 

under the assumption that z k (t } s;() are known. The double time labels in A (t, s,h; 
and a r (t,s;C) distinguish these symbols as dynamical objects, in contrast to the non- 
dynamical quantities A (t,h; and a (t,() bearing single time arguments in (3.30). In 
particular, (3.31) shows that semiclassical admissibility is preserved under evolution. 
Using methods already employed in §3a we directly compute 



A(t, s) = A w (f , h; - V H )r(s, t)e">-' = A w (*, h; - V J exp {tu • Z(t, s)} 

l u = I 

= exp{Z(t,s).V}AJt,h;0). (3.32) 



u = 
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Taking the Wigner transform, using the cluster expansion theorem and isolating the clas- 
sical flow yields 

oo _. 

AJt, s, h; C) = e ^w(M,^;0- v exp j ^ _^ (( . ^ ^ ^ y) j ^ ^ n . ^ a | C) ) . (3.33) 

In the above formulae V is understood to act on the symbol A w (t } U; • ). The cluster oper- 
ator Cj contains j such gradients. Hence if A w (t, U; () is polynomial in ( the ^2^=2 ^ jl ^ 
series will reduce to a finite sum, and the expansion of e SZw '^ will similarly terminate. 
This kind of simplification underscores the fact that most of the nontrivial dynamical 
effects are built into Z w (t } s, U; (). 

Formula (3.33) verifies the Ti-anarytic structure of the evolving symbol A (t,s) which 
was asserted in (3.31). Moreover, by extracting the coefficients of common powers of 
h in (3.33) one obtains explicit formulae for a in terms of the known quantities 

{z (t (t, s; () , a k (t, (*)}£ = o • The first three of these are found to be 

a Q (t,s;() = A c (t,g(t,s\()), (3.34) 

a x (t, s; C) = a x (t, g(t,s\()) + ^(t, s 5 • VA C (<, g(t,s\()) , (3.35) 



,(t, s; C) = a 2 (t, g(t,s\C)) + 2z 1 (t } s; () • Va, (t, g(t,s\()) + (z^t, s; Q-Vf 



(3.36) 



+z 2 (t, s; C) • V - ^B 12 2 X g(t,s\() • V + ^B 12 B 23 X g(t,s\() • v] A c (t, g(t,s\()) . 



Clearly this method of determining a r requires no further integration of ODE's beyond that 
needed to determine the z r coefficients. In this sense {z r } act as a 'basis' with which to 
represent arbitrary evolving observables. 

(n) 

Representations of the ^-coefficients 7^ / of the Heisenberg-Weyl evolution map T s t 
are obtained from (3.33) by choosing the observable symbol to be a static h independent 
phase space function, say a. Then a (t, • ) = 8 a and 



7™{a) = a(g{t,3\-))=g{t,sy 



while from (3.35) 



7$(a) = z i(M; •) ' Va(flf(*,s| •))• 



(3.37) 
(3.38) 



For a Weyl quantized system, (3.36) leads to 



7 $(a) = z 2 (t,s; ■)-Va(g(t,s\-)) -^B 12 2 ^g a (t,s), gP(t,syV a Vpa(g(t,s\-)) 

+ ^B 12 B 23 ^g a (t } s) , /(*, s) , g"<(t, s)yV a V p V ^a(g{t,s\ ■ )) . (3.39) 
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Note that, by the linearity of T g t , in order to compute A (t,s) it is really sufficient to 
determine just {7^0}. For then (3.31) follows from (3.30) and (1.8) with 

ar(*^;-) = EQ^ )a r-„(*>-)- 

This approach will be taken in the next section. 

The derivation leading to (3.33) is formally exact — no approximations have been 
made nor terms neglected. The result provides a representation of evolution which explic- 
itly exposes the underlying classical motion and systematically characterizes higher order 
quantum corrections in terms of a cluster-graph expansion. Representation (3.33) is uni- 
versal in the sense that it holds for all Hamiltonians H(t) and observables A(t) that are 
semiclassically admissible operators. Identity (3.34) states that to leading order in h the 
evolving Weyl symbol A w (t } s } U; () is given by the classical observable pulled along the 
classical flow. This is the analog, for Weyl symbols, of Egorov 's TheoremJ 42 ] In (3.33) 
the dynamical quantum corrections to classical behavior are consolidated into two aspects: 
1) SZ w which accounts for the deviation of the quantum trajectory Z w about its classical 
limit g, and 2) the h dependence of the cluster operators C- which has its origin in the 
*-product. 



IV — Moyal Bracket Method 

Here a pure classical transport recurrence relation method of computing the semiclas- 
sical expansion (3.31) is constructed. This alternative approach does not rely on the Z w 
representation (3.3) or the cluster expansion of §2b. The transport based methods of this 
section are an extension of techniques found in the prior literature^ 9-11 ' 39 ' 43 ] on semiclas- 
sical expansions rooted in the Moyal formalism. Typically these works have assumed that 
the Hamiltonian is static. For the more general case of explicitly time dependent operators 
which we wish to consider, certain simplifying features of the static case are absent. In 
these circumstances it is found that the backward evolution equation is the appropriate 
Heisenberg equation of motion. Subsection 4a explains this and summarizes the properties 
of evolution operator T t . In the next subsection the alternative method is developed and 
consistency with the results of §3 is considered. 

4a. Heisenberg and Moyal Equations 

The essential features of the various possible Heisenberg equations of motion which 
need to be clarified are all present at the operator level, so we set this discussion there. 
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Let A(t), and H(t), be a time dependent Schrodinger picture operator, and Hamiltonian. 
In the Heisenberg picture, 



A(t, s) = r(s, t)A(t) = U(t, s)l A(t) U(t, s), (4.1) 

and H(t,s) = r(s,t)H(t). Note that if H(t) is not static, then U(t,s) is not e (*-*)#M 
and it need not commute with H(t), so H(t,s) ^ H(t). The 'forward' and 'backward' 
Heisenberg equations of motion are, respectively, 

d - 1 - dA 

di A(t ' s) = m s) ' H(t ' s) ] + li {t ' s) ' (4 - 2a) 

^-A(t,s) = -^-[A(t,s),H(s)}, (4.26) 

where ^(t,s) = r(s,t)(dA(t)/ dt). These follow from the four variants of the Schrodinger 
equation (1.3) in which either U or is differentiated with respect to either its first or 
second argument (in deriving these, ?7(t,s)^ = ?7(t,s) _1 = U(s,t) is used). 

The following features of the Heisenberg equations are critical. The backward equation 
(4.26) involves the commutator with the Schrodinger picture Hamiltonian H(s). It is clearly 
a homogeneous linear non-autonomous evolution equation for A(t, s) with s playing the 
role of the evolution parameter. Time t is passive, and forms a convenient 'pinning time' 
for the initial condition 

A(t,T)\ r=t = A(t). (4.3) 

By contrast the forward equation is not a homogeneous evolution equation for A(t, s): 
the Schrodinger operator dA(t)/dt is an arbitrary prescribed function used to form an 
inhomogeneous term. Moreover, the commutation operator - , H(t,s)\ involves 

the Heisenberg picture Hamiltonian. Therefore both coefficients H(t,s) and ^(t,s) are 
unknown (since U(t, s) is), and their dependence on the initial time s violates the standard 
formulation of an evolution equation. A proper evolution equation is, however, obtained 
from (4.2a) in the special case where A and H are static. 

Thus, it is the backward Heisenberg equation which provides a solvable statement 
of the evolution problem in the general case of time dependent operators. An intuitive 
understanding of this fact is also useful. Equation (4.1) shows that r(s,t) maps an object 
'at' time t, namely A(t), into one at time s (recall that A(t,s) is traced with states frozen 
at time s to yield expectation values). In this sense then, the initial time parameter of the 
Schrodinger equation is in fact the natural forward evolution parameter for the Heisenberg 
picture. We now denote it by r, as was done in (1.1). For other applications of the 
backward Heisenberg equations see Refs. [44,45]. 
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The Wigner-Weyl image of the preferred Heisenberg equation (4.26) is the Moyal 
equation (1.1). Introducing T Tt as in (1.7) it reads 

^ T /a) = -{r T /a),HjT)} M , (4.4a) 
for arbitrary phase space functions a. The companion initial condition (4.3) becomes 

T M (a) = a. (4.46) 

(identity (4.4a) appeared in the symbolic form (1.12).) 

Note that (4.4a) contains the Moyal bracket rather than the Poisson bracket. Roughly, 
this means it involves infinite order differentiation of T Tt (a) [cf. (2.18)] and so it is 
not a first order transport equation like (1.2). A rigorous analysis^ 23 ] would interpret 
{ • , if w (r)} M as a pseudodifferential operator acting in a Banach space of functions, say 
T 3 a. Under suitable conditions^ 46 ! on this operator, (4.4) will be a linear evolution 
problem with a unique solution defined for r £ R. We shall proceed under the assumption 
that this is the case. A parallel situation exists regarding the operator valued Schrodinger 
equation (1.3). Restrictive assumptions on H(t) must be imposed in order to ensure the 
applicability of a theory of linear evolution equations in 7i in a form that guarantees the 
existence of a unitary solution. For example, if the self-adjoint family of operators H(t) 
has a stable domain and on this domain is strongly continuously differentiable in r then 
this suffices to prove (e.g. Refs. [47,48]) the existence of a unique U(t } s). 

Clearly, 1^ t is a kind of fundamental solution for the Moyal equation. As the next 
lemma shows, it possesses all the properties characteristic of an evolution family. 



Lemma 4. (Properties of T s t ) For all s,t £ R 

(a) T g t : T ^> T is C-linear. 

(b) T g t is bijective and T s = T t . 

(c) An evolution law is obeyed: for all p £ R 

r s , ( = r v or^ ( . (4.5) 

(d) r M = idp 

(e) If 1 £ T is the constant function with value 1, then r f l = 1. 

(f) r t preserves the *-product: for all a, 6 £ T 

r s , t (a)*r Sjt (b) = r Sjt (a*b). (4.6) 
Hence T s t preserves the Moyal bracket too. 

If T g t is defined by (1.7) then the proofs of these properties follow routinely from 
properties of U(t,s). Alternatively, one can define T g t (a) as the unique solution of (4.4). 
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The proof of Lemma 4 then relies in a central way on this unique solvability hypothesis, 
and is not difficult so we omit it. The index order chosen for T s t is in harmony with 
evolution law (4.5). Taken together, properties (a), (b) and (f) state that each T s t is an 
automorphism of the algebra T with product *. 

The above properties of T s t do not display any feature which is parallel to the unitarity 
of the Schrodinger evolution U(t,s). It is, however, possible to formulate the problem in 
such a way that T s t is unitary. Let B 2 denote the Schmidt class of operators^ 49 ! acting in 
U = L 2 (R d ). It is knownt 50 ! that B 2 is a Hilbert space with respect to the inner product 
(A, B) i — ^ Tr(A^ B). From the unitarity of U(t, s) it is simple to check that the Heisenberg 
evolution r(s,t) — acting in B 2 — is unitary. Now, the Wigner transform a is a Hilbert 
space isomorphism from B 2 onto T = C 2 (R 2d , h~ d d() [cf. (5.2a)]. Thus, it follows that 
the Heisenberg- Weyl evolution map t = a w o r(s,t) o cr" 1 is unitary in this space T of 
square integrable symbols. 

An interesting application of the properties above is to verify preservation of the 
canonical Moyal brackets (1.10). 

Corollary. Let Z^(t,s) = T st ( y 7r z ^, where ir z are Euclidean coordinates on R 2d , i.e. 
7r* = c w (5*) is projection onto the z th component. Then (1.10) holds. 

Proof: Since 7r* supports only one nonzero derivative, {tt 1 , vr J } M = {tt 1 , 7r J } = J tJ . 
From this and properties (a), (e) and (f) one computes 

{zi(t,s), zi(t,s)} M = r M ({A tt'-} m ) = Tjri) = r*. 

■ 

An advantage of the Weyl symbol representation is that it is possible to identify 
a natural ^-parity of quantum mechanics. Commonly occurring physical systems have 
Weyl symbols (3.1) that are finite order in h with all h-o&d terms absent. In this stan- 
dard situation H w (t } U; () has even ^-parity and, as the lemma below shows, so does the 
Heisenberg- Weyl evolution operator T s t . This simplifies the semiclassical expansion (1.8) 
since all 7^ with odd n will vanish. 

Lemma 5. (Evolution ^-Parity) Assume the Hamiltonian's Weyl symbol H w (t } U; •) 
has even dependence on fi. Then Y s t is even in ft,, and so all jg 2 ^ +1 ^ = 0. 

Proof: The map t (where the ft dependence has been made explicit) is defined as 
the unique solution to (4.4). If ft is replaced by —ft throughout these equations, it is found 
that r~^* also obeys (4.4) because { • , • } M is ft-even (cf. (2.18)) and H w (t) is assumed to 
be so. By uniqueness of the solution, t = . ■ 
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4h. Classical Inhomogeneous Transport Equations 

Consider the Heisenberg-Moyal evolution problem (4.4). The Hamiltonian symbol 
H (t) is given by (3.1). We take (1.8) as an ansatz for the analytic h structure of the 
Heisenberg-Weyl evolution map Y s t . The initial conditions on 7^ induced from (4.46) 
are 

7 t ( ?=Id^, 7* ( ? = (n>l). (4.7) 

If these two representations are substituted into (4.4a) and the Moyal bracket h expansion 
(2.18) is used, then an ^-expanded version of the Moyal equation results. Equating in it 
the coefficients of h r , r > 0, leads to the following identities which 

a r (M; - )=7#(a) (4.8) 

must satisfy: 

a 2 a r (t, r; () + {a r (<, r; • ) , H c {t, •)}(() = f r (t, r, (), (4.9a) 

j+fc + 2/ = r ^ / " •/ " 

< r 

Here d 2 is the partial derivative on the second time argument, and the sum is over all 
integers j } k } l > obeying the two indicated constraints. By definition f = 0. Equation 
(4.9a) is an inhomogeneous linear Hamiltonian transport equation. For r = it reduces to 
the homogeneous equation (1.2). The inhomogeneous term / is built from a Q , . . . , a r _ 1 , 
which leads to recursive solutions for 7^ . The unique solvability of (4.9) is based upon 
the following results. 

Lemma 6. (Backward Flow) The flow g(t,s\-) of a time dependent Hamiltonian 
H c (t, ■ ) obeys 

d 2 g(t,r\0 = -Vg(t,T\OJVH c (r,0- (4.10) 



Proof: Differentiating the flow identity 

g(T,t\g(t,T\())=( (4.11) 

with respect to r using the chain rule and solving for d 2 g yields 

d 2 g(t,r\() = -V^t,^)- 1 ^^,^), 

where z = g(t,T\(). Now differentiate (4.11) with respect to ( to find V(/(r, t|z) _1 = 
V(/(t,r|C). Substituting Wg from here and d 1 g from Hamilton's equation (1.6) results in 
(4.10). Compare also Lemma 4.2.32 of Ref. [13]. ■ 
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Theorem 2. Let the Hamiltonian H c (t } () } inhomogeneous term /(r, (*) and initial 
data /3(C) be smooth functions, with g the flow of H . Then the unique smooth solution 
6 of the Hamiltonian transport problem 

db(r, C) + {6(r, • ) , H c (t, •)}(() = /(r, C) (4.12a) 

b(t,C) = m (4-126) 

is given by 

&(r,C) = /%(t,r|C))- / ^/(p,^,r|C)). (4.13) 



Proof: [Uniqueness] Let b be a solution of (4.12). Upon substituting ( = g(r, t\z) into 
(4.12a) it follows from Hamilton's equation that the left side assumes the form of a total 
derivative, so (d/dr) 6(r, g(r, t\z)^ = f(r, g(r, t\z)^ . Integrating this ODE from r to i and 
using (4.126) yields 

b(r,g(T,t\z)) = P(z) - / dpf(p,g(p,t\z)). 



Finally, choose z = g(t } T\() to find that b must be given by (4.13). 

[Existence] If b is defined by (4.13) it clearly obeys (4.126), so it must be shown to satisfy 
(4.12a). Differentiating (4.13) with respect to r and using (4.10) to represent d 2 g leads to 

<96(r,C)= [-V/%(t,T|C))-V^,r|C) 

+ j* dp Vf(p, g(p,r\0) • Vg(p, r |0] JVH c (r, () + f(r, (). (4.14) 

By the chain rule, V • Vg = dfl(g(t, t|(*)) /d(, with a similar reduction for V/ • Vg in 
the integral term. The first term on the right of (4.14) thus reduces to a Poisson bracket, 
specifically — {6(r, •), H c (t } •)}. Thus 6 obeys (4.12a). ■ 

Application of Theorem 2 to the transport problem (4.9) and (4.7), of quantum origin, 
immediately yields the solutions 



a r (t,T;() = 6 rtQ a(g(t,T\())+ £ 



(") 



^ 2 2l (2i + i)\j\k\ 

j+k + 21 = r v ' 

k<r (4.15) 



x J dpB 2l+1 ^R k (t,p; •), hjip, ■)y{g(p,T\0). 

For r = the integral term is absent and a Q (t,r;C) here agrees with (3.37). For r > 1, 
equations (4.15) provide a recursive integral solution based on classical transport. 
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It is instructive to see the explicit form of the first few coefficients a r . We consider 
a system with a Weyl quantized Hamiltonian; accordingly h ■ = 0. The parity property of 
T g t assures us that all the odd index coefficients & n vanish. The predictions of (4.15) for 
the non- vanishing terms to order U 6 are: 

a 2 (M;C)= / dp — BUa (t,p; -),H c (p, -)y(g(p,s\()) 

•J S 

a 4 (M;C) = dp^BUa (t,p; - ),H c (p, ■)y-^B 3 ^ 2 (t,p; ■ ), H c (p, ■ )^| (g(p,s\()) 
a 6 (M;C) = J dp^^BUa (t,p; - ),H c (p, ■ )y + ^5 5 ^a 2 (t, p; - ),H c (p, -)y 

-^B 3 ^(t,p; -),H c (p, ■)yyg(p,s\()). (4.16) 

Further simplifications of the recursive representation (4.15) arise if H c (t } () is a poly- 
nomial function of order N. In this case H can support at most N non- vanishing deriva- 
tives, so all factors with B 2l+1 -<a k (t, /?;•), H (p, ■ )>- are zero if 2/ + 1 > N. An important 
special case of this arises if H c (t } () is quadratic. As is evident from (4.16), & n = 
for n > 0, and the zeroth order semiclassical expansion (1.8) is then exact. Specifically, 

r = 7 (0 ? 

L s,t ls,t ■ 

The iterated Poisson bracket transport formulae (4.15) and (4.16) are remarkably 
concise representations. This simplicity is somewhat deceptive due to their recursive na- 
ture. The computation of a 2 in terms of the original symbol a, for example, requires the 
substitution of a Q = a o g into the first of (4.16) with subsequent evaluation of its third 
derivative. Further iteration to compute a 4 , etc., will lead to more derivatives and multiple 
time integrals 

(n) 

Representations of the Heisenberg-Weyl semiclassical expansion coefficients 7^ / are 
obtained from (4.15) by making the notational change indicated in (4.8). For example, the 
first correction is 

7m«(C)= / dp{a{g(t,p\-)),h 1 (p, ■)}{g(p,s\()) (4.17a) 

J s 

= {a(g(t,s\-)), J dph 1 {p,g(p,s\-))}(Q, (4.176) 

where invariance of the Poisson bracket under the symplectic map g was used. Similarly, 
equations (4.16) and (4.8) immediately provide recursive forms of 7^ , j^t anc ^ l^t m 
the case of a Weyl quantized Hamiltonian. 

We close this subsection by verifying that representation (4.17) of 7^ agrees with the 
previous result (3.38). To do this, simply substitute (3.28) into (3.38) and note that the 
derivative structure V •a((/)V J (/* = V -(a o g), again by the chain rule. In this way (4.176) 
emerges. 
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V — Semiclassical Behavior of Expectation Values 



The preceding sections have developed, in considerable detail, methods for computing 
the small h expansion of Heisenberg picture operators. From a physical standpoint it is 
ultimately desirable to employ these results in the broader context of observable dynamical 
phenomena. In particular, a semiclassical treatment of quantum expectation values is 
called for. Viewed in this way, the h expansion of A (t,s) forms one individual theme, 
which should be combined with an appropriate description of quantum states to yield a 
semiclassical analysis of the expectation value. In this section we consider these questions 
in a general way. Our objective here is not to construct detailed higher order expansion 
formulae, but rather it is twofold: (1) to sketch in rough outline the nature and advantages 
of the theory of expectation values in the Heisenberg picture and Weyl representation; (2) 
to indicate the evident physical character of the results which can be obtained. 

In general a quantum state is a density operator p } that is, a positive trace-class 
operator of unit trace, Tr p = 1. The expectation value of observable A in the state p is 
(A), = Tr(Ap). The Weyl symbol versions of these concepts are well known. The trace 
becomes a phase space integral 

Tr A = h~ d f d(A w (Ch (5-1) 

as follows at once from (2.6). To obtain the trace of a product one can combine Groe- 
newold's formula (2.9) with the identity 

j d(f(V 1 .JV 2 )-<A w ,B w y(0 = j d( f(0)AJOBJO 

(which follows for analytic / by formally integrating V 2 by parts); thus one finds 

Tr (A(t, s)p) = h~ d J d( AJt, s; OpJO- (5.2a) 

We summarize a few elementary analytical properties^ 2 ' 3 ' 51-54 ] of densities. The symbol 
p is a real valued, uniformly continuous L 2 function on R 2d , which vanishes at infinity 
and is bounded by \p w \ < (2/h) d . The integral h~ d J d( p (C) = 1 corresponds to the unit 
trace normalization, via (5.1). Unlike a classical probability density, the quantum density 
p can take on negative values. 

An important class of states are the pure states. These are 1-dimensional projection 
operators. If tp £ 7i is a normalized vector in the range of such a pure state /3, then 
p = IV'KV'I an d (5.2a) becomes 

<^(M))^ = j d(A w (t,s;()w^(), (5.26) 
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where = h~ d \tp) (ip\ w is the 'Wigner function' of tp. 

Of critical importance to semiclassical analysis is an understanding of the h analytic 
structure of pure state densities. The projection property of pure states, p 2 = p, in symbol 
form reads p * p = p or using Groenewold's formula, 

e W^^p w , Pw y = Pw . (5.3) 

Suppose, for example, that one may represent the structure of p near h = by a Laurent 
series, 

oo 
j = M 

where either M £ Z (the integers) or M = — oo, and substitute this series into (5.3) and 
match common powers of h. It will be found that either M = or M = — oo. The only 
regular (M = 0) series consistent with (5.3) is p =1, which is not a density (it is not 
normalizable). Thus M = — oo and we conclude that the Wigner function of a pure state 
(or a finite linear combination of pure states) cannot be regular at h = 0, and if it has 
an isolated singularity there, it must be an essential singularity. Note that this conclusion 
holds regardless of whether or not the wave function ij)(x) representing the pure state has 
explicit h dependence. 

By contrast, mixed states (p 2 ^ p) can have Ti-regular densities. An example is 
the density p{fi) = e~@ H /Tr e~@ H of the canonical ensemble at inverse temperature fi. 
In this case the Ti-regular expansion of p(f3) is the venerable Wigner-Kirkwood^ 2 ' 55-58 ] 
expansion. This expansion can be derived by proceeding as in §2d after replacing t — > fih/i. 

Note that equations (5.2) have been written in the Heisenberg picture, where the time 
evolution is ascribed to the observable A(t,s) and the state, p or ?/>, is held fixed. This 
is in keeping with the expansions (3.33) and (4.15). Much of the prior literature^ 39 ' 46 ' 59 ] 
uses the Schrodinger picture in which the density (initially p) carries the time dependence 
according to p(t,s) = U(t, s) p U(t, s)^ = r(t,s)p. Here the roles of t and s are opposite 
to the observable case (1.4), and p(t,s) obeys the quantum Liouville equation of motion 

^p(t,s) = ^[H(t),p(t,s)] (5.4) 

(usually considered for static Hamiltonians where H(t) = H). 

Clearly, one could try to apply our results for the Heisenberg- Weyl evolution map 
to compute p (t,s) = T t s (p w )- However, this Schrodinger picture approach encounters 
considerable inherent difficulties due to the fact, noted above, that typically the initial 
density p is essentially singular in h. Clearly one would have to determine the precise 
nature of the essential singularity as it evolves in time, and develop techniques to construct 
the small h perturbative form of the residual nonsingular factor, if one is to get useful 
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semiclassical expansion formulae. This is one of the reasons we are led to the Heisenberg 
picture, where the observable maintains a regular (semiclassically admissible) h analytic 
form. 

In the circumstances where the density p carries an h essential singularity, the deter- 
mination of the semiclassical form of the expectation value (A(t,s))^ requires a careful 
evaluation of the phase space integral (5.26). We investigate this topic, in the preferred 
Heisenberg picture, by considering several informative examples. 

Example 1. (Semiclassical States and Interference) Consider a state vector ip of semi- 
classical form in the configuration representation, 

ip(x) = <p(x)exp \—S(x)\, (5.5) 
n 

with smooth amplitude <p> £ L 2 (R d ) and smooth real phase function S. Combining (5.26) 
and (2.7), the Wigner function based expectation of observable A(t,s) in this state is 



{A(t,s)) = dq(27rn)- d dpdvA w (t } s } n;q } p)p(q + v/2)p(q-v/2) 

J J J (5.6) 

X exp { % - [ - p • v + S(q + v/2) - S(q - v/2)] } . 

For each fixed q } we will find the leading h J, asymptotics of the (p, u)-integral in (5.6) 
by employing the stationary phase asymptotic formula [Ref. [60], Theorem 2.2] 

{2TT7i)- n ' 2 j d n yF(y } n)e^/ % ^^ (5.7) 
= |detVV$(y o )|" 1/ V(y o ,^)e( J / ?t )*(^) e J ( 7r / 4 ) si s nVV *(^) + 0{n 1 ' 2 ). 

In this formula y o denotes the critical point (V$(y ) = 0) of the smooth phase function $, 
whose Hessian is VV$. In applying (5.7) to (5.6), identify y = (p, v), hence n = 2d } and 
$(p, v) = —p • v + S{q + v/2) — S{q — v/2), while F(p } v } h) is the non-exponential portion 
of the integrand. The unique critical point of $ is easily found to be given by 

P = VS(q), v o = 0, (5.8) 

while the Hessian there is 



w$( Po ,o 



6 
6 



(5.9) 



It is readily demonstrated that this matrix has the eigenvalues +1 and —1, each of multi- 
plicity d. Hence the signature is zero. Concerning A w (t } s, U; q,p), we will keep its leading 
semiclassical value a Q (t, s; q } p) as given by (3.34). Employing (5.7)-(5.9) in this way gives 

<i(M))^= / dqA c (t,g(t,s\q,VS(q)))mq)\ 2 +0(n 1 / 2 ). (5.10) 
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Formula (5.10) has a transparent physical interpretation, insofar as IVK?)! 2 is the 
probability density at q and VS'(g) is the leading order momentum of tp. Note that while 
the Wigner function has an essential singularity in h, nevertheless once the momentum 
integration is completed the result is an ^-regular expectation value. 

It is of interest to generalize this example by considering a state which is a (normalized) 
superposition of two semiclassical wave packets like (5.5), 

The corresponding density 

P = Pi + h + IV'iX^ I + I^XV'i I 

is still a pure state, (?/>|, but the cross terms will now describe interference of the two 
packets. In computing ^ A(t, s) )^ the terms p x and p 2 will each give a contribution of the 
form (5.10). The stationary phase analysis in (5.6)-(5.10) readily extends to the cross term 
integrals. The rapidly oscillating phase is now $(p, v) = —p • v + S^q + v/2) — S 2 {q — v/2), 
with critical point v o = and p o = |V(S' 1 + S 2 )(q)- Again the Hessian VV$(p o ,v ) has 
determinant ( — ) d and signature 0. To within 0(h 1 ^ 2 ) the cross term contribution is found 
to be 

< i ^ 5 )> cr oss = / ^A c (t, 5 (t,3|g,|V(5 1 + 5 2 )(g)))2Re^ 1 ^ e ^-^)/ ?t ]( (Z ). (5.11) 

The surviving essential singularity, with a phase difference S 1 — S 2 , is characteristic of 
interference phenomena. Formula (5.11) is suitable for computing the effects of wave 
packet interference at the physical value of h. It should be kept in mind that <p 1 and 
(p 2 are initial, not evolving, amplitudes. Of course, if h J, this term oscillates to zero 
whenever S 1 7^ S 2 . 

The computation of higher order h corrections to the leading terms, (5.10) and (5.11), 
would require a more detailed asymptotic analysis of the Wigner function and the phase 
space integral. This is beyond our scope here, which is to indicate the character of the 
obtainable result. A more comprehensive investigation of the h structure of physically 
relevant quantum densities p w is a significant topic in its own right, and is the subject of 
ongoing work! 38 ' 61 ' 62 ] in the literature. The fact that in (5.10) and (5.11) the formal phase 
space symmetry between q and p is broken is merely the result of employing a configuration 
space wave function ij)(x) in these examples. 

Example 2. (Harmonic Oscillator) Consider the 1-dimensional harmonic oscillator 
with unit mass, H = (p 2 + uj 2 q 2 )/2, and choose the initial state to be the h independent 
Gaussian 

<f>(x) = ce" x2 / 2 , c = tt- 1 / 4 . (5.12) 
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For illustrative purposes we pick the observable A = q 2 , and shall compare the computation 
of (q 2 )^ , ip t = U(t,0)(f), via two methods: (i) standard propagator based quantum 
mechanics, and (ii) the Heisenberg-Weyl formalism. Because of the simple quadratic 
Hamiltonian, both of these approaches yield closed form solutions. 

(i) The lowest order WKB approximation^ 63-68 ! gives exact representations for the prop- 
agators of (q } p)-quadratic Hamiltonians. For the harmonic oscillator the resulting expres- 
sion is known as Mehler's formula^ 69 ' 70 ] 



K(t,x,y) = (x\U(t,0)\y) 



uj/h 



sin cut 



e -z(7r/4V(i) exp | ^2 + y 2^ CQt ut _ 2xy cgc ut 



(5.13a) 

valid for t ^ mr/co (n £ Z), with the Maslov index /i being constant between successive 
conjugate points: /i(t) = 2n + 1 for mr < cut < (n + 1)tt. At a caustic the propagator is a 
distribution, 

K( — ,x,y)=e- in *' 2 6((-) n x-y). (5.136) 

cu 

The fact that /i(t) varies discontinuously in going through a caustic is a necessary 
feature of the propagator, which ensures that K is the kernel of the t-strongly continuous 
evolution operator U(t, 0). The fact that /i depends only on t is a special feature of static 
quadratic Hamiltonians, such as the harmonic oscillator. 

Calculating for t ^ mr/uj } the evolved state is obtained as a complex Gaussian integral 



A 2 esc 2 cut 

i\ cot cut 



A 



sin cut 



\ 2 

e -i(^/4)M*) (i _ iX C ot cut) " 1/2 exp { X — 



1 — iX cot cut 



where A = to /Ti > and the branch cut for ( • ) -1 / 2 is along the negative real axis with 
argument in (— 7r,7r). A careful check will show that i^ t {x) has a removable singularity at 
cut = mr and is indeed continuous there. 

Once tp t is used to calculate an expectation value, (^4.)^ = J dx tp t (x)[A'ip t ](x), any 

effect of the Maslov index will be lost. We take the example A = q 2 and find after a 
Gaussian integration that 

<? 2 >* t =\{^ sm 2 cut + cos 2 cut). (5.14) 

Initially this result holds away from caustics, but a separate calculation for t = mr/co using 
(5.136) establishes (5.14) for all t £ R. Note that the final result is regular in h and t. 

(ii) For initial state <f> of (5.12) the Wigner function is simply computed as the Fourier 
transform of a real Gaussian, 

„2 



9 h Ju ™ 
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and possesses an essential singularity at h = 0. This must be paired with the evolved 
symbol A w (t,h; q, p). As was shown in §4b, the quadratic nature of H implies that T t 

reduces to its classical part 7^ = g(t \ ■ )*. Computation of the classical harmonic oscillator 
flow g = (q c ,p ) is trivial; one finds q c (t\q } p) = (p / to) sin cot + qcoscot. The expectation 
value is then 

<9(^0) 2 )^= / / dqdp q c (t\q, pf w^{q,p) = -(—sin 2 tot + cos 2 tot), 

in agreement with (5.14). 

The examples above, and others^ 43 ! from the literature, point to three viable compu- 
tational methods for extracting the semiclassical behavior of an observable's expectation 
value, i) Calculate (if possible) the phase space integral (5.2) exactly. This was the ap- 
proach used in Example 2. ii) Use the stationary phase technique applied to the (p, v) 
integration in (5.6), in order to obtain the leading small Ti form of <( A(t, s) } . This was 
done in Example 1. iii) Employ a hybrid evaluation of the expectation value. Specifically, 
in (5.2) keep the exact form of the initial density matrix p w (() intact with h fixed at its 
physical value, but replace A w (t } s; () with its semiclassical expansion (3.31), truncated at 
some finite order in fi — most often after the leading term A c (t, g(t, s|C)) . This approach 
has the merit of being easy to implement in practical applications. At the same time it 
respects those aspects of the uncertainty principle (at the physical value of Ti) which are 
encoded into the initial quantum state. 

In a sense, the third approach appears to treat the h dependence in A (t, s; dif- 
ferently from that in /? w (0? but it is a consistent procedure since A (t, s; is Ti-regular, 
while p (C) is typically a localized uniformly bounded L 1 (R 2d ) function. The mathemati- 
cally distinct asymptotics problem implicit in the hybrid method may at first sight appear 
unnatural, but it does have a sound physical basis. The limit conventionally described 
by L 7i I 0' is physically achieved when some ratio U n /S n becomes small, where S is some 
quantity characteristic of the system. Since h is fixed at its physical value, S must become 
large, e.g. by varying the observable or Hamiltonian. By doing this while leaving the initial 
state fixed, one realizes the asymptotics inherent in the hybrid approach. 

Further insight into the Moyal semiclassical method results from a comparison with the 
WKB approximation, I 35 ' 70-72 ! which provides a semiclassical description of the propagator 
K. The propagator argument Q = (t,s,x,y) defines a set of distinct classical trajectories 
{z-{ ■ ; Q)} obeying qj(s) = y and qj(t) = x. Depending on Q } the solution set to this two 
point boundary value problem may be empty, finite or infinite. Locally the trajectories 
depend continuously on Q } but globally Q-space is fragmented by surfaces where the 
nature of the solutions changes. The WKB propagator ansatz^ 70 ~ 72 ^ is a 'sum over paths' 
of contributions from each trajectory z ■ of the form 

K j(Q) = (27rn)~ d/2 \D j (Q)\ 1/2 e~ zin/4)tlliQ) exp{^Sj(Q) + 0(fi 1 )}, 
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where S-, D- and // • are, respectively, the principal function (action), its VanVleck de- 
terminant and a Maslov-related phase for the j th trajectory. The solvability analysis of 
the Hamiltonian two point boundary value problem and the determination of the Maslov 
index changes which accumulate through caustics along a trajectory (where D ■ becomes 
singular) are, for general Hamiltonians, significant mathematical problems. Their solution 
easily becomes prohibitively difficult for practical calculations. 

In contrast to this situation, the Heisenberg-Weyl semiclassical method is based upon 
the classical flow, i.e. the unique classical trajectories which exist for each initial value 
( £ R 2d . Furthermore, there is no need to compute Maslov indices for these trajectories. 
In summary, it would seem that for many purposes, while the ( x | U(t, s) | y ) based compu- 
tations are correct, the use of this intermediate tool can introduce significant unnecessary 
complications due to caustics. 

VI — Conclusion 

The Heisenberg picture evolution operator r(s,t), and equivalently its Weyl phase 
space counter part T s t , constitute a complete description of quantum evolution. The 
Heisenberg-Weyl operator T s t — mapping symbols of observables into symbols — admits a 
regular small h expansion that provides a fundamental statement of semiclassical dynamics. 

Explicit forms of the operator valued expansion coefficients 7^ are found by a two 
stage process. First, the dynamical phase space coordinate operators Z(t,s) = r(s,t)z, 
via their symbol Z (t, s, U; (*), are obtained by using cluster-graph (Theorem 1) and Green 
function (Lemma 3) methods to formally solve their equation of motion (1.9). The h 
expansion of Z (t, s,h; () has as its leading term the exact classical flow g(t } s\() with 
higher order h terms accounting for quantum fluctuations about the classical motion. 
Having computed the quantum trajectory Z (t, s, U; () the next stage is to determine the 
evolution of an arbitrary time dependent observable. Again an application of the cluster 
expansion yields an exact representation, (3.33), of the dynamically varying observable. 
Parsing this representation according to powers of h constructs the semiclassical expansion 
operators 7^. 

As the explicit expressions (3.37)-(3.39) for 7^ show, the differential structure of 
these operators is a universal function of the Poisson bracket operators and finite or- 
der phase space gradients. The effect of a specific choice of Hamiltonian H is confined to 
the classical flow arguments g(t } s\() and the value of the higher order trajectory fluctua- 
tions z r (t,s;(). For all fi-even Hamiltonians the odd coefficients jg 2 ^ +1 ^ vanish, thereby 
increasing the accuracy and efficiency of the expansion. For (*-quadratic Hamiltonians the 
classical approximation is exact, namely T s t = 7^°^. 
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Prior literature on the semiclassical expansion (1.8) of the Heisenberg-Weyl evolution 
operator is evidently sparse. However, for the expansion of the Weyl symbol of a time 
evolving observable, as in equation (3.31), there are a variety of results. The conceptual 
ingredients (found in Section 4) of employing transport recurrence relations obtained from 
an order-by-order h expansion of the Moyal equation of motion goes back to Berezin and 
ShubinJ 9 ' 11 ] When applied to densities rather than observables, the recurrence transport 
method has been treated extensively by Prosser J 10 ' 39 ] 

The most widely studied semiclassical expansion in the literature is the WKB approx- 
imation for the propagator (x | U(t,s) | y ). As such, it establishes a natural benchmark 
for measuring the effectiveness of the Heisenberg-Weyl expansion (1.8). In the expansion 
(1.8) there is only one classical trajectory determined by an initial value problem. The co- 
efficients 7^ are singularity free and, if the classical motion does not run away to infinity 
in finite time, are defined for all time. Unlike the WKB approximation there is no essential 
singularity in h, no need to solve a two-point variational problem for the relevant classical 
trajectories, and no small time restrictions or caustic difficulties. 

It was emphasized in §3 that this globally unique trajectory aspect of expansion 
(1.8) is the joint result of employing Weyl symbols for Heisenberg picture evolution. It 
is worthwhile to have a somewhat heuristic understanding of this theme. Certainly the 
use of phase space based symbols (such as Weyl's) admits the possibility of a unique 
trajectory being the relevant one (in the Weyl case its initial value is (*, the symbol variable). 
Such a possibility seems naturally excluded, e.g., for the usual (x \ ■ ■ ■ \y) representation 
because the pair of configuration variables (x } y) is not equivalent to a point in phase 
space. However, use of the Weyl symbol alone in studying quantum evolution is not 
sufficient to guarantee absence of caustics. For example, if the Schrodinger picture is used 
and U w (t,s) is studied, the relevant boundary condition^ 38 ' 40 ] for trajectories is {z c (t) + 
z c (s)) /2 = (*, which admits caustics and carries no global guarantee of unique solutions. 
Why is the situation so favorably different for the Heisenberg picture evolution operator 
r(s,t)t Roughly, this seems to be because, through the appearance of both and U, 
the Heisenberg operator r(s,t)A maps between states at a common time, s, which turns 
out to be the initial time for the trajectory value (. 



APPENDIX A: AFFINE CANONICAL COVARIANCE 



The Wigner-Weyl picture of quantum mechanics possesses a simple covariance with 
respect to affine canonical transformations.! 11 ' 16 ' 20 ' 26 ' 73 ' 74 ^ The first part of this appendix 
reviews the nature of this covariance, while the second part examines precisely how this 
property is respected by the formulae and methods of computation derived in §3. 
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Let $ : R 2d — > R 2d be an affine canonical transformation. For convenience we represent 

$ as 

*{z)=S~\z-z \ (Ala) 



or equivalently 



$- 1 (0 = S( + z o , (A.lb) 



where the symplectic matrix S £ Sp(2d) and the displacement z o e U 2d are arbitrary. Let 
us say that a unitary operator V corresponds to $ if the following quantum counterpart 
of (A. la) holds, 

V z 1/ t = $(5) = S _1 {z - zj). (A.2a) 

In this case 

V* zV = Sz + zJ. (A.2b) 

Given $, we note that such corresponding operators V = V(S } z ) always exist. Indeed 
they may be represented by V = T(z o )M(S) where T(z o ) is the Heisenberg phase space 
translation operator^ 16 ' 75 ' 76 ] 

iX^exp^.J-^), (A3) 

and M(S) is a metaplectic operator J 16 ' 26 > 77 ] These operators respectively obey the funda- 
mental identities T(z )^ z T(z ) = z + z I and M(Sy zM(S) = Sz. Condition (A. 2a) may 
be shown to determine V uniquely up to a phase (i.e. if V , V obey (A. 2a) then V = cV 
with |c| = 1). 

In terms of such corresponding classical and quantum affine transformations, the co- 
variance of the Wigner transform takes the form: 

(V*AV) w = A w o $-i (A4) 

for any operator A Before proving this identity, it is of interest to clarify its interpretation. 

Regard $ as a (passive) change of phase space coordinates, zh(e &( z )- If A w is a 
classical observable, then A w o $ _1 represents the same observable in the new coordinates. 
On the other hand, the unitary operator V can be used to implement a change of basis 
vectors in 7{, defined by tp \— > cf> = V^tp. Under this 'quantum coordinate change' the 
observable A transforms into V^AV. Thus (A. 4) says that, for corresponding transforma- 
tions, the coordinate change process commutes with the Wigner- Weyl mapping between 
quantum and classical observables. 

We now give the proof of (A. 4). Conjugating the representation (2.4) of A in terms 
of its symbol, with V , results in 



Vl 'AV = (2vr)- 2d J du J d(A w (()e tu - 



iu ■ (Sz+z — C) 
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Upon making the change of integration variables v = S T u and w = S 1 (( — z ) one finds 

V^AV = (27r)- 2d J dv j dwAJSw + z o )e lv - ( "- w \ 
which implies (A. 4). 

Let us now consider the covariance of the semiclassical method developed in §3. For 
brevity, denote conjugation by a unitary V with a superscript, A v = AV . From a gen- 
eral perspective, the transformation properties of all semiclassical symbols are immediately 
given by (A. 4); thus 

(i(t,3)^) w (C) = A w (t,3,^;$- 1 (C)). (A5) 
Taking the example A = z and combining this with the h expansion (3.3) yields 

(Z(t,sr)JO = Y / - T ^s^- 1 (0). (A.6) 



r 

r=0 



Assume that the affine canonical map $, corresponding to V , is h free (i.e. that S and z o 
don't depend on U) } as is natural from a semiclassical perspective. Then (A. 6) will be a 
bona fide h expansion of the form 



C(M,^;C) = ]L^z^(M;C), ( A - 7a ) 

r=0 



oo ^ 



with 

v 



(t,3;C) = z r (t,3;$- 1 (C)) (A.7b) 



being the covariance relation for the coefficient functions. In the remainder of this ap- 



pendix, we indicate how property (A. 5) emerges if A(t,s) is computed by the methods 
of §3, assuming that (A. 7) is given. In this regard, we remark that it is possible to derive 
(A. 76), independently of the above discussion, by examining the Jacobi equations obeyed 
by in an inductive fashion. We will not present this inductive proof since it is tedious 
and several of its facets would be duplicated by the proof of (A. 5) which we sketch next. 
To begin the derivation, conjugate the operator equation (3.32) with V , 

A(t, s) v = exp {Z(t, s) v • V}A w (t, %■ 0). (A8) 

Note that it is A w [not (A v ) w ] which appears on the right side, because operators [not 
symbols] are being conjugated. Using the cluster theorem, the Wigner transform of (A. 8) 
is 

oo 1 

Al(t, s, h; C) = e^w(M^;0- v exp j £ £ v (C; ^ ^ ^ y) j ^ ^ % . Z ^ (M | C) ) , (A9 ) 
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in which the relevant cluster operators are, cf. (3.86), 

Cj((; h, t, s, V) = ]T J] [ e («/2)B„ _ !] x Z^(*, 3, n ; C) • V. (A.10) 

CeCj aEE(C) 3 

Consider a typical extended bracket B = V fc • JV^ acting on the & th and / th instances 
of Z^ in the product function. Notice that (A. 7) implies at once the version of (A. 5) 
restricted to the case A = z, viz. 

z^(t } s } n ; o = z w (t } s } n^- 1 (o). 

Using this identity to replace each Z^ in (A. 10) leads, by the chain rule, to the replacement 
of B with 

a 

Vfc-VS-^OJVS-^CfV,, (All) 

where now V fc acts on the original symbol Z (t, s, h; ■ ), followed by evaluation at $ -1 ((*). 
Since $ _1 is canonical (V$ _1 = S is symplectic) the bracket (A. 11) reduces to precisely 
B . Moreover, since $ _1 is affine, the derivative S is ( independent and so the above 
bracket invariance extends to arbitrary higher order derivatives arising from any e ( tfl / 2 ) B a 
in (A. 10). This results in the covariance rule for the cluster operators, 

Cj((; h, t, s, V) = C 3 (*-\Q; h, t, s, V) . (A12) 

Finally, use of (A. 7b) and (A. 12) to replace the V dependent objects on the right side 
of (A. 9) completes the verification of (A. 5). The 7i coefficient covariance ajf = a r o $ _1 
follows from this. 



APPENDIX B: INTEGRAL KERNEL REPRESENTATIONS 

Like most quantum operators, the Heisenberg-Weyl evolution r(s,t) has an integral 
kernel realization. In this appendix we present two established phase space representations 
of Heisenberg-Weyl evolution and clarify their linkage to T s t . 

The kernel representation of T s t is 

r * *(«)(*)= / d(K(s,z\t,()a(0- O^ 1 ) 

Formally, K is the distribution valued solution of (1.12) which has the s = t delta function 
initial value £> z (()- The ordering of the time labels in (BA) reflects the fact that s is the 
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active time variable in the backward equation of motion (4.4a) and that t is the initial 
pinning time. Given the existence of K, setting a(() = S ,(() shows that 

T a>t (6 z ,)(z) = K(3,z\t,z'). (B.2) 

Formula (B.2) provides the link between the T s t operator and the kernel K. Putting 
a = 8 Z into (3.37)-(3.38) explicitly displays the distributional nature of the formal h J, 
asymptotic expansion of K. 

One standard form of K is realized in terms of A(z) } the quantized delta function 
operator^ 78 ' 79 ] defined by A(z) = h d S z . The work of Grossmann^ 80 ] and Royer^ 54 ! has 
shown that A is essentially the product of two elementary unitary operators, namely 
A(z) = 2 d T(2z)H. Here T is the Heisenberg translation operator (A. 3), and II is the 
parity operator, Hz = — zll. The usefulness of A is apparent from the identities 

A = J Y d a ^ A ^' < z ^ = T 4 Aa ( z )} ■ (B.3a, b) 

The first of this pair provides another way to Weyl quantize the symbol a, comparable 
to (2.4) or (2.5). The second identity determines the Weyl symbol a from the operator 
A via the trace, and is an alternative to the Wigner transform (2.7). Straightforward 
demonstrations of (-B.3) may be found in Ref. [53]. 

We now review an exact representation^ 36 ' 53 ] of K(s, z\t, (). Employing (B. 3a) allows 
one to write the Heisenberg evolution of A as 

T(M)i = / ^(M)tA(CMM)A w (C)- (BA) 

The symbol of r(s,t)A is T g t (A w )(z). A second form of this symbol is obtained from the 
right side of (BA) by multiplying by A(z) and taking the trace to give 

K(s,z\t,() = h- d Tr[U(t,syA(()U(t,s)A(z)]. (B.5) 

This formula states that kernel K(s, z\t, () is the expectation value that results from the 
s—*t evolution of observable A(()/h d with respect to the 'density matrix' A(z). The 
kernel K is real valued, and satisfies the symmetry condition K(s, z\t, () = K(t, (\s } z), as 
is easy to show from (B.5). It also obeys a composition law implied by (4.5). 

As the recent work of Marmot 81 ] has revealed, K admits a path integral representa- 
tion. For a Weyl quantized Hamiltonian, the standard time slicing method in combination 
with the short time approximation U w (t + At,r)(z) ~ exp ( — iAtH c (r, z)/h) leads to 

K(s,z\tX) = j- d f DzDieWW'd. 
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Here Dz and are infinite products of the dimensionless measures h d dz and h d d^. 
The phase is the action-like functional 



dr 



i(r) • j'z{r) + Jf c (r,S(r) + ±f(r)) - H c (t,z(t) - ±f(r)) 



The phase space paths z( ■ ) obey the boundary conditions z(s) = z, z(t) = (*, whereas the 
variations £( • ) are unrestricted. 
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Abstract 

The Moyal description of quantum mechanics, based on the Wigner-Weyl isomorphism 
between operators and symbols, provides a comprehensive phase space representation of 
dynamics. The Weyl symbol image of the Heisenberg picture evolution operator is regular 
in h and so presents a preferred foundation for semiclassical analysis. Its semiclassical 
expansion 'coefficients,' acting on symbols that represent observables, are simple, globally 
defined (phase space) differential operators constructed in terms of the classical flow. The 
first of two presented methods introduces a cluster-graph expansion for the symbol of an 
exponentiated operator, which extends Groenewold's formula for the Weyl product of two 
symbols and has h as its natural small parameter. This Poisson bracket based cluster 
expansion determines the Jacobi equations for the semiclassical expansion of 'quantum 
trajectories.' Their Green function solutions construct the regular h J, asymptotic series 
for the Heisenberg- Weyl evolution map. The second method directly substitutes such a 
series into the Moyal equation of motion and determines the h coefficients recursively. In 
contrast to the WKB approximation for propagators, the Heisenberg- Weyl description of 
evolution involves no essential singularity in h, no Hamilton-Jacobi equation to solve for 
the action, and no multiple trajectories, caustics or Maslov indices. 



